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Massive neutrinos open up the possibility for a variety of new phys¬ 
ical phenomena. Among them are oscillations and double beta de¬ 
cay. Furthermore they influence several fields from particle physics 
to cosmology. In this article the concept of massive neutrinos is 
given and the present state of experimental research is extensively 
reviewed. This includes astrophysical studies of solar, supernova 
and very high energy neutrinos. Future perspectives are also out¬ 
lined. 
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1 Introduction 


The birth of the neutrino due to W. Pauli in 1930 was a rather desperate 
attempt to explain the continuous /3-spectrum [1]: 

” I have considered ... a way out for saving the law of conservation of energy. 
Namely, the possibility that there could exist in the nuclei electrically neutral 
particles, that I will call neutrons (which are today called neutrinos) which 
have spin 1/2 and follow the exclusion principle. The continous /-spectrum 
would then be understandable assuming that in /-decay together with the elec¬ 
tron, in all cases, also a neutron is emitted in such a way that the sum of 
energy of neutron and of electron remains constant... I admit that my solution 
appears to you not very probable... But only who dares wins, and the gravity 
of the situation in regard to continuous /-spectrum...” 

The experimental discovery of the neutrino by Cowan and Reines [2] in 1956 
and the observation that there exist different types of neutrinos by Danby 
et al. [3] were important milestones. The last important step about neutri¬ 
nos stems from the LEP-experiments measuring the Z°-width which results 
in 2.993 ± 0.011 flavours for neutrino masses below 45 GeV [4]. 
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From all particles of the standard model, neutrinos are the most unknown. 
Because they are treated as massless particles, the physical phenomena as¬ 
sociated with them are rather limited. On the other hand in case of massive 
neutrinos , which are predicted by several Grand Unihed Theories, several new 
effects can occur. This article reviews the effects of massive neutrinos as well 
as the present knowledge and experimental status of neutrino mass searches. 


2 Theoretical models of neutrinos 


The presently very successful standard model of particle physics contains 
fermions as left-handed chiral projections in doublets and right-handed 
charged fermions as singlets under SU(3)cG)SU(2)£,®U(l)y transformations. 
Neutrinos only show up in the doublets which does not allow any Yukawa 
coupling and therefore no mass with the minimal particle content of the stan¬ 
dard model. Moreover, because neutrinos are the only uncharged fundamental 
fermions, they might be their own antiparticles. 

In the following chapter, a theoretical description of neutrinos is given as well 
as possible extensions of the standard model to generate neutrino masses. A 
second requirement will be to explain why neutrino masses are so much smaller 
than the corresponding charged fermion masses. The most promising way is 
given by the see-saw-mechanism. 


2.1 Weyl-, Majorana- and Dirac-neutrinos 


The neutrino states observed in weak interactions are neutrinos with helicity 
-1 and antineutrinos with helicity -|-1. For massless neutrinos and the absence 
of right-handed currents there is no chance to distinguish between Dirac- and 
Majorana neutrinos . Because V-A theory is maximal parity violating the 
other two states (neutrinos with helicity -|-1 and antineutrinos with helicity 
-1), if they exist, are unobservable. If neutrinos are massless a 2-component 
spinor (Weyl-spinor) is sufficient for description, hrst discussed for the general 
case of massless spin 1/2 particles by Weyl [5], which are the helicity -1(-|-1) 
projections for particles (antiparticles) out of a 4-component spinor T. They 
are given by 



(1 ± 75)4/ 


( 1 ) 


The eigenvalues of 75 (chirality) agree with those of helicity in the massless 
case. Here the Dirac equation decouples into two seperate equations for 
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respectively. An alternative 2-coniponent description was developed by Majo- 
rana [ 6 ] to describe a particle identical to its antiparticle. If neutrinos acquire a 
mass, in general both helicity states for neutrinos and antineutrinos can exist, 
making a 4-component description necessary. Here a 4-component Dirac-spinor 
can be treated as a sum of two 2-component Weyl-spinors or as composed out 
of two degenerated Majorana neutrinos. However it is still an open question 
whether neutrinos are Dirac or Majorana particles. The Majorana condition, 
for a particle to be its own antiparticle, can be written as 

C-\C = -7j (2) 

with C as the charge conjugation operator. The real charge conjugated state 
ji)^ is not obtained by a C operation but by CP, because pure charge 
conjugation results in the wrong helicity state. In the case of a Dirac-neutrino, 
the helds 'ipR and are sterile with respect to weak interactions and therefore 
they are sometimes called Nr and . The most general mass term in the 
Lagrangian including both Dirac- and Majorana helds is given by 

^ + ^l^r) + "if + rriR'ipL'^R) + h.c. 


with 
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In the general case of Ua active neutrinos and n* sterile neutrinos Ad is a 
{ua + Us) X (na + Us) matrix (see [7]). Assuming only one neutrino generation, 
diagonalisation of Ai results in the eigenvalues 

"ii,2 = + "if) =t {m¥ — rnuY + (4) 


Four different cases can be considered: 

- mf = mf = 0 —> mi ^2 = "i^: The result is a pure Dirac-neutrino which 
can be seen as two degenerated Majorana helds. 

- , mf —> mi ^2 ~ "i^: Neutrinos are Pseudo-Dirac-Neutrinos. 

- = 0 — mi ^2 = rn¥ ,mR'- Neutrinos are pure Majorana particles. 

- thr 3> m^,m¥ = 0: This leads to the see-saw-mechanism . 

The see-saw-mechanism [8,9] results in two eigenvalues 
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Because neutrino masses should be embedded in GUT-theories, the latter 
offers two scales for and m^. All fermions out of a multiplet get their 
Dirac-mass via the coupling to the same Higgs vacuum expectation value. 
Therefore the neutrino Dirac mass is expected to be of the order of the charged 
lepton and quark masses. The heavy Majorana mass can take values up to the 
GUT-scale, which is in the simplest models about 10^® GeV. Assuming three 
families and a unique the classical quadratic see-saw 




: = ml : ml : 


(7) 


emerges. This is only a rough estimate because several effects influence this 
relation. Instead of the quark-masses, the charged lepton masses could be 
used. In scenarios where is proportional to m^ for the different families, a 
linear see-saw relation results. Depending on the GUT-model, the mass scale 
of m^ need not be related to the GUT-scale but might be in connection 
with some intermediate symmetry breaking scale (Table 1). Last not least 
the relation holds at the GUT scale, to get a prediction at the electroweak 
scale, the evolution has to be calculated with the help of the renormalisation 
group equations. Especially the third term can experience a signihcant change 
depending on the used GUT model like [11] 

= O.OSm^ : 0.09m^ : 0.38m^ SUSY-GUT (8) 
■ rnu^ = O.OSm^ : 0.07m^ : 0.18m^ SO(10)-GUT (9) 

A further see-saw-mechanism resulting in almost degenerated neutrinos is dis¬ 
cussed in chapter 2.3. 


2.2 Massive neutrinos in the standard model 


In the present standard model with minimal particle content, neutrinos remain 
massless. The simplest extension to create neutrino masses is the inclusion of 
SU(2) singlet states denoted by N^. Because of hypercharge zero they remain 
singlets of the entire gauge group and have no new interaction with the gauge 
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Table 1 


Predictions for neutrino masses according to different see-saw models. As can be 
seen, the quantity {ruu^) , measured in neutrinoless double beta decay , corresponds 
in most cases to rriu^ (after [ 10 ]). 


model 

my. 

{my,) 

my^ 

my^ 

Dirac 

1-10 MeV 

0 

100 MeV-1 GeV 

1-100 GeV 

pure Majorana (Higgs triplet) 

arbitrary 

my. 

arbitrary 

arbitrary 

GUT seesaw (M lO^^^GeV) 

10-1^ eV 

my. 

10-® eV 

10-® eV 

Intermed. seesaw (M ~ 10®GeV) 

10“^ eV 

my. 

10-2 eV 

10 eV 

SU(2)®SU2®U(1) (M 1 TeV) 

10-1 eV 

my. 

10 keV 

1 MeV 

light seesaw (M <C 1 GeV) 

1-10 MeV 

<C rUy, 



charged Higgs 

< leV 

<C rUy, 




bosons. New Yukawa-couplings of the form 

C = K{iyLeL)(^_']NR +h.c. (10) 


result in a Dirac mass term of = huV 2 where V 2 ~ 246 GeV is the vacuum 
expectation value of the neutral component of the standard model Higgs- 
doublet. In order to produce an eV-neutrino , the Yukawa-coupling hy has to 
be smaller than 10“^'^. Some hne-tuning is required for this, on the other hand 
the generation of the mass pattern is still unknown and such a small hy might 
be possible. An immediate consequence of a mass term is, that similar to the 
quark sector, a mixing between the mass eigenstates z/j and flavour eigenstates 
z/q, can occur 




( 11 ) 


allowing several new phenomena, e.g. neutrino oscillations , which will be 
discussed later. Nevertheless the global lepton number L remains a conserved 
quantity. 

Without introducing additional fermion singlets, it is only possible to generate 
Majorana mass terms, because only ur and its charge conjugate (z/l)*" exist. 
These terms necessarily violate L and therefore also B-L by two units. The 
only fermionic bilinears carrying a B-L net quantum number are 

iiCT)? . («)Ar (12) 

The necessary extensions of the Higgs-sector to produce gauge invariant 
Yukawa couplings therefore offer three possibilities: a) a triplet b) a single 
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charged singlet and c) a double charged singlet. 

Case a: The additional Higgs triplet A = (A°, A“, A ) carries hypercharge 
-2 and the neutral component develops a vacuum expectation value of v^. It is 
this vacuum expectation value which enters the Yukawa-coupling for the mass 
generation of neutrinos . There is no prediction for the masses or V 3 , but it can 
be much smaller than V 2 and therefore explain the lightness of neutrinos. This 
additional vacuum expectation value would also modify the relation between 
the gauge boson masses to [12] 



mlcosWw 1 + 44 
'^2 


(13) 


which, by using experimental values, results in 

— < 0.07 (14) 

V 2 


Case b: This corresponds to the Zee-model [13]. By introducing a single 
charged higgs /i_ and additional higgs doublets, Majorana masses can be gene¬ 
rated at the one-loop level by self-energy diagrams. If only one higgs couples 
to leptons, a mass matrix of the following form can be derived [14] 


M = mo 


0 a cosa 
a 0 sina 
cosa sina 0 


(15) 


with 


tana = A(1 - 


fe 


mt 


a = 


e/i ml 


cosa 


er 
2 . 


mo = Ami fer / cosa 


(16) 

( 17 ) 

(18) 


where / are the Yukawa coupling constants and the electron mass is neglected. 
This in general implies two nearly degenerated neutrinos and one which is 
much lighter. 

Case c: By including an additional double charged higgs /;;++ with (B-L) quan¬ 
tum number 2, it is possible to generate masses on the 2-loop level which are 
therefore small [15]. It can be shown that for three flavours one eigenvalue is 
zero or at least much smaller than the others. 
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All the solutions described above violate B-L by introducing B-L breaking 
terms in C. On the other hand, the vacuum could be non-invariant under B- 
L, for example as a spontaneous breaking of a global B-L symmetry. This is 
discussed in more detail in connection with the associated Goldstone boson, 
called majoron, in chapter 4.1. 


2.3 Neutrino masses in grand unified theories 


As already seen in the description of the see-saw-mechanism, by choosing 
a large it is possible to get small neutrino masses. To hnd a scale for 
, an implementation of this mechanism into grand unihed theories seems 
reasonable. The simplest grand unihed theory is SU(5) even if the minimal 
SU(5)-model is ruled out by proton-decay experiments. Because all the fun¬ 
damental fermions can be arranged in one multiplet there is no room for a 
right-handed neutrino and consequently no Dirac-masses. Minimal SU(5) is 
also B-L conserving which is given by the multiplets and the gange invariance 
of the higgs held conplings. For this reason Majorana mass terms also do not 
exist. Therefore in the minimal SU(5) nentrinos remain massless. By extend¬ 
ing the higgs-sector it is possible to create mass terms via radiative corrections 
as in the Zee model. Nevertheless the proton decay bound remains. 

The next higher grand nnihed theory relies on SO(IO). All fnndamental 
fermions can be arranged in a 16-mnltiplet, where the 16th element can be 
associated to a right-handed nentrino . This allows the generation of Dirac 
masses. In SO(10) B-L is not necessarily conserved opening the chance for 
Majorana mass terms as well. The breaking of SO(IO) allows diherent schemes 
like 


AO(IO) ^ SU{5) SU{3) ® SU{2) ® U{1) (19) 


or into a left-right symmetric version after the Pati-Salam model [16] 


SO{10) ^ SU{2)l® SU{2)r® SU{4) (20) 


This generates a right-handed weak interaction with right-handed gange 
bosons. These models create neutrino mass matrices like [17,18] 


^ fvL \ 

fvR j 


and vl 


A(m^)L 

vr 


( 21 ) 


where f is a 3x3 matrix and vl,vr are the vacnnm expectation valne of the 
left-handed and right-handed higgses respectively. Diagonalisation leads to 
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masses for the light neutrinos of the form 




fKmw)L Dr-l 


vr 


-m f {m ) /vr + 


( 22 ) 


Two important things emerge from this. First of all, the first term domi¬ 
nates over the second, the latter is corresponding to the quadratic see-saw- 
mechanism . Because no neutrino masses are involved in the hrst term and if 
f is diagonal, no scaling is included resulting in a model with almost degen¬ 
erated neutrinos in leading order. This is sometimes called type II see-saw- 
mechanism [17]. In case the hrst term vanishes, we end up with the normal 
see-saw-mechanism . For a more extensive discussion on neutrino mass gener¬ 
ation in GUTs see [12]. 


3 Kinematical tests of neutrino masses 


3.1 Beta decay 


The classical way to determine the mass of Fe is the investigation of the electron 
spectrum in beta decay. A hnite neutrino mass will reduce the phase space 
and leads to a change of the shape of the electron spectra, which for small 
masses can be investigated best near the Q-value of the transition. In case 
several mass eigenstates contribute, the total electron spectrum is given by a 
superposition of the individual contributions 

N{E) a F{E. Z)-p-E-(Q-E)-J2 1(Q - Ef - m? | Ul \ (23) 

2=1 


where F(E,Z) is the Fermi-function, the m* are the mass eigenvalues and 
are the mixing matrix elements. The different involved m* produce kinks in 
the Kurie-plot where the size of the kinks is a measure for the corresponding 
mixing angle. This was discussed in connection with the now ruled out 17 keV 
- neutrino [19,20]. 

The search for an eV-neutrino near the endpoint region is complicated due to 
several effects [21,22]. The number of electrons in an energy interval AE near 
the Q value scales with 


n{Q — AE) oc ( 

8 


AE 

~Q 


3 


(24) 
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Fig. 1. Mainz 1998 raw data prior to publication (private communication by Ch. 
Weinheimer). The signal/background ratio is increased by a factor of 10 in com¬ 
parison with the 1994 data. Eq corresponds to the center of mass of the rota¬ 
tion-vibration excitations of the molecular ground state of the daughter ion ^HeT~^. 

making a small Q-value advantageous, but even for tritium with the relatively 
low endpoint energy of about 18.6 keV only a fraction of 10“® of all electrons 
lies in a region of 20 eV below the endpoint. A further advantage of tritium 
is Z=l, making the distortions of the (3 - spectrum due to Coulomb - interac¬ 
tions small and allow a sufficiently accurate quantum mechanical treatment. 
Furthermore, the half-life is relatively short and the involved matrix element 
is energy independent (the decay is a superallowed transition between mirror 
nuclei). All these arguments make tritium the favoured isotope for investiga¬ 
tion. 

For a precise measurement, the resolution function of the used spectrometer 
has to be known quite accurately. Additionally also the energy loss of electrons 
in the used source, consisting of molecular tritium T 2 , is important. Effects 
of molecular binding have to be taken into account and only about 58 % of 
the decays near the endpoint lead to the ground state of the ^He T’''-ion, 
making a detailed treatment of hnal states necessary. A compilation of the 
obtained limits within the last years is given in Table 2. As can be seen, most 
experiments end up with negative £t values, which need not to have a 
common explanation. For a detailed discussion of the experiments see [21,22], 
While until 1990 mostly magnetic spectrometers were used for the measure¬ 
ments, the new experiments in Mainz and Troitzk use electrostatic retarding 
spectrometers [23,24], Fig. 1 shows the present electron spectrum near the 
endpoint as obtained with the Mainz spectrometer. The negative values 
for a larger interval below the endpoint are understood for both experiments. 
While in the Troitzk experiment, using a gaseous T 2 source, the energy loss 
of trapped electrons in the spectrometer was underestimated, for the Mainz 
experiment, using a thin him of T 2 , roughening transitions in the him seem to 
be the reason. More recently, the Troitzk experiment observed excess counts 
in the region of interest, which can be attributed to a monoenergetic line short 
below the endpoint. This is currently under study in the Mainz experiment 
which after some upgrades might explore a Ve mass region down to 2 eV. 
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Fig. 2. spectrum obtained with a cryogenic bolometer. The big spikes corre¬ 

spond to calibration peaks (from [26]). 

A complementary result would be the measurement of /3-decay in . Be¬ 

cause of its endpoint energy of only 2.6 keV, according to eq.(24) it allows 
a high statistics search near the endpoint. A cryogenic bolometer in form of 
a Re-foil together with a NTD-germanium thermistor readout has been suc¬ 
cessfully constructed and a measurement of the /3-spectrum above 100 eV was 
obtained (Fig. 2) [25]. Because this experiment measures the total released 
energy reduced by the neutrino rest mass, energy loss and final state effects 
are not important. 

CPT-invariance assures that rriy^ = rriy^ . A direct measurement of rriy^ as 
proposed by [27] is the internal bremsstrahlungs - spectrum in EC-processes 

Z -\- e —> (yZ — 1) -1- + T (25) 


The most convenient isotope is and the limit obtained is [28] 

< 225eR {95%CL) (26) 

This is rather weak in comparison with beta decay. Astrophysical limits on 
will be discussed in chapter 7. 


3.2 Pion decay 


The easiest way to obtain limits on m,y^ is given by the two-body decay of the 
7 i~^. For pion decay at rest the neutrino mass is determined by 

= ml+ + ml+ - 2m^+{pl+ + m2+)(F2) (27) 

Therefore a precise measurement of the muon momentum and knowledge 
of and is required. These measurements were done at the PSI resulting 
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Table 2 

Compilation of existing /3-decay results of tritium and the corresponding Vf. mass 
limit. _ 


experiment 


mp, (eV) 

Tokyo (INS) 

-65 ± 85 ± 65 

< 13.1 

Los Alamos (LANL) 

-147 ±68 ±41 

< 9.3 

Zurich 

-24 ±48 ±61 

< 11.7 

Livermore (LLNL) 

-130 ± 20 ± 15 

< 7.0 

Mainz 

-22 ± 17 ± 14 

< 5.6 

Troitzk 

1.5 ±5.9 ±3.6 

< 3.9 


in a limit of [29] 

= (-0.016 ± 0.023)Mel^=^ ^ < 170A;e\^(90%C'L) (28) 

where the largest uncertainty comes from the pion mass. Investigations of pi- 
onic atoms lead to two values of = 139.56782 ± 0.00037 MeV and = 
139.56995±0.00035 MeV respectively [30], but a recent independent measure¬ 
ment supports the higher value by measuring = 139.57071 ±0.00053 MeV 
[31]. 

3.3 Tau-decays 


Before discussing the mass of Vr h should be mentioned that the direct de¬ 
tection of Vr via CC reactions still has not been observed and all evidences 
are indirect. The goal of E872 (DONUT) at Fermilab is to detect exactly 
this reaction. With their presently accumulated data of 4.55 ■ 10^'^ protons 
on target, about 60 Vr CC events should be observed. The present knowl¬ 
edge of the mass of Vr stems from measurements with ARGUS (DORIS II) 
[32], CLEO(CESR)[33], OPAL [34], DELPHI [35] and ALEPH [36] (LEP). 
Practically all experiments use the r-decay into five charged pions 

r —> ± 57r^(7r°) (29) 


with a branching ratio of BR = (9.7 ± 0.7) ■ 10“^. To increase the statistics 
CLEO, OPAL, DELPHI and ALEPH extended their search by including the 
3 TT decay mode. But even with the disfavoured statistics, the 5 prong-decay is 
much more sensitive, because the mass of the hadronic system peaks at about 
1.6 GeV, while the 3-prong system is dominated by the Oi resonance at 1.23 
GeV. While ARGUS obtained their limit by investigating the invariant mass 
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Fig. 3. 2-dimensional plot of the hadronic energy versus the invariant mass of the 
5(6) vr-system. The error ellipses are positively correlated, because both the hadronic 
mass and the hadronic energy are determined from the momenta of the particles 
composing the hadronic system (from [36]). 

Table 3 

Comparision of Ur mass limits (95 % CL) as measured by various experiments. 
Numbers with * include also events from 3 7r-decay. 


experiment 

number of events 

rriy^ limit (MeV) 

ARGUS 

20 

31 

CLEO 

266 

30 

OPAL 

2514* -h 5 

27.6 

DELPHI 

6534* 

27 

ALEPH 

2939* -h 41 

18.2 


of the 5 TT-system, ALEPH, CLEO and OPAL performed a two-dimensional 
analysis by including the energy of the hadronic system (Fig. 3). A hnite 
neutrino mass leads to a distortion of the edge of the triangle. A compilation 
of the resulting limits is given in Table 3 with the most stringent one given by 
ALEPH [36] 


< 18.2MeH(95%C'L) (30) 

Plans for a future charm-factory and B-factories might allow to explore rriy^ 
down to 1-5 MeV. 

Independent bounds on a possible Vr mass in the MeV region arise from pri¬ 
mordial nucleosynthesis in the early universe. Basically, three effects influence 
the detailed predictions of the abundance of light elements [37]. An unstable 
Vr or its daughters would contribute to the energy density and therefore influ- 
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ence the Hubble-expansion. Moreover, if they decay radiatively or into 
pairs, they would lower the baryon/photon ratio. A decay into hnal states 
containing v^. or z/g would influence the neutron fraction and therefore the 
abundance. Recent analysis seems to rule out Dirac masses larger than 
0.3 MeV and Majorana masses larger than 0.95 MeV at 95 % CL for [38]. 
An independent constraint from double beta decay , only valid for Majorana 
neutrinos, is discussed in chapter 4.1. 


4 Experimental tests of the neutrino character 


Double beta decay 


The most promising way to distinguish between Dirac and Majorana neutrinos 
is neutrinoless double beta decay . For extensive reviews see [39-41]. Double 
beta decay was first discussed by Goeppert-Mayer [42] as a process of second 
order Fermi theory given by 

(Z, A) —(Z + 2, A) + 2e“ + 2z7g decay) (31) 

and subsequently in the form of 

{Z, A) —> (Z + 2, A) + 2e“ {0u(3(3 decay) (32) 


hrst discussed by Furry [43]. Clearly, the second process violates lepton number 
conservation by 2 units and is only possible if neutrinos are massive Majorana 
particles as discussed later. In principle V+A currents could also mediate 
neutrinoless double beta decay , but in gauge theories both are connected 
and a positive signal would prove a hnite Majorana mass [44,45]. To observe 
double beta decay, single beta decay has to be forbidden energetically or at 
least strongly suppressed by large angular momentum differences between the 
initial and hnal state like in ‘^^Ca . Because of nuclear pairing energies, all 
possible double beta emitters are gg-nuclei and the transition is dominated 
by O’*" —O'*" ground-state transitions. The 2z//5/5 decay can be seen as two 
subsequent Gamow-Teller transitions allowing only virtual l’''-states in the 
intermediate nucleus, because isospin selection rules forbid or at least strongly 
suppress any Fermi-transitions. The matrix elements for the 2vj3j3 decay can 
be written in the form [41] 


., 2 . ^(0;i|t-a||l+)(l+||Aa||0+) 

Y E, + Q/2 + me-E, 


(33) 
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Fig. 4. Different spectral shapes of observable sum energy spectra of emitted elec¬ 
trons in double beta decay . The n=l,3,7 forms correspond to different majoron 
accompanied modes, n=5 (grey) is the 2i/(5(5 decay and the Qufifi decay results in a 
peak. The energy values are taken for . 

and for the Oz//3/3 decay as 


^GT = ■ crn||0+) 

m^n 

(34) 


(35) 


with t_ as the isospin ladder operator converting a neutron into a proton, a 
as spin operator, r =\rm — \ and H{r) the neutrino potential. In the neu¬ 

trinoless case because of the neutrino potential also other intermediate states 
than 1+ might be populated [46]. 

Typical energies for double beta decay are in the region of a few MeV dis¬ 
tributed among the four leptons which are therefore emitted as s-waves. The 
phase space depends on the available Q-value of the decay as oc (Oz//3/3 de¬ 
cay) and oc (2z//3/3 decay), numerical values can be found in [47]. From the 
experimental point of view, the sum energy spectrum of the two emitted elec¬ 
trons has a continuous spectrum for the 2h'j3j3 decay , while the Oz//3/3 decay 
mode results in a peak at the position corresponding to the Q-value of the 
involved transition (Fig.4). The single electron spectrum for the two nucleon- 
mechanism is given by [48] 

(Ti, T 2 , COS0) = (Ti + 1)2(T2 + 1)2 X 6 {Q - Ti - T 2 )(l - pip 2 Cose) (0z//3/?)(36) 

(Ti, T 2 , cosO) = (Ti + 1)2(T2 + 1)2 X (g - Ti - T2)^(1 - (3iP2Cose) (2z//3/3)(37) 

where Ti,T 2 are the kinetic energies in units of the electron mass, j3i^2 is the 
velocity and 6 the angle between the two electrons. Some favourite isotopes 
are given in Table 4. 
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Experimental considerations A rough estimate of the expected half-lives 
for the 2vj3j3 decay mode results in the order of 10^° years. Therefore it is an 
extremely rare process making low-level counting techniques necessary. To ob¬ 
tain reasonable chances for detection, isotopes with large phase space factors 
(high Q-value) and large matrix elements should be used. Also it should be 
possible to use a signihcant amount of source material, which is improved by 
second generation double beta decay experiments using isotopical enriched ma¬ 
terials. One of the main concerns is possible background. Background sources 
are normally cosmic ray muons, man-made activities like , the natural 

decay chains of U and Th, cosmogenic produced unstable isotopes within the 
detector components, and . The cosmic ray muons can be shielded 

by going underground, the natural decay chains of U and Th are reduced by 
intensive selection of only very clean materials used in the different detector 
components, which is also valid for , and by using a lead shield. To avoid 
cosmogonies, the exposure of all detector components to cosmic rays should 
be minimized. This is important for semiconductor devices. can be re¬ 

duced by working in an air-free environment, which can be done by using pure 
nitrogen. For more details on low-level connting techniqnes see [49]. 

The experiments focusing on electron detection can be either active or passive. 
Active detectors have the advantage that source and detector are identical as 
in the case of , but only measure the sum energy of both electrons. On 
the other hand passive detectors allow more information like measuring en¬ 
ergy and tracks of both electrons seperately, but usually have smaller source 
strength. Under the assumption of a flat backgronnd in the peak region, the 
sensitivity for the Qvj3j3 half-life limit can be estimated from experimental 
quantities to be 


TO; 


1/2 


oc a\ 


I M-t 

B-AE 


(38) 


where a is the isotopical abundance, M the used mass, t the measuring time, 
AE the energy resolution at the peak position and B the background index 
normally given in counts/year/kg/keV. Some experiments will be described 
in a little more detail. 

Semiconductor experiments : In this type of experiment , hrst done by Fiorini 
et al. [50], source and detector are the same, the isotope under investigation 
is . The big advantage is the excellent energy resolntion (typically abont 
5 keV at 2 MeV). However, the techniqne only allows the measurement of 
the snm energy of the two electrons. A big step forward was done by using 
enriched germaninm (natnral abundance of : 7.8 %). The Heidelberg- 
Moscow experiment [51] in the Gran Sasso Laboratory is using 11 kg of Ge 
enriched to 86 % in form of hve HP-detectors. A background as low as 0.2 
connts/year/kg/keV at the peak position has been achieved. To improve fur¬ 
ther on background reduction, a pulse shape analysis system was developed 
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to distinguish between single site events (like double beta decay) and multiple 
site events (like multiple Compton scattering) which seems to improve B by 
another factor of hve. The IGEX collaboration is using about 6 kg in form of 
enriched detectors [52], 

Moreover, there is always the possibility to deposit a double beta decay emit¬ 
ter near a semiconductor detector to study the decay, but then only transitions 
to excited states can be observed by detecting the corresponding gamma rays. 
Scintillator experiments : Some double beta decay isotopes can be used as part 
of scintillators. Experiments were done with in Form of CaF 2 [53] and 
in Form of CdW 04 [54]. 

Cryogenic detectors'. A technique which might become more important in the 
future can be bolometers running at very low temperature. Such detectors 
normally have a very good energy resolution. At present only one such exper¬ 
iment is running as a 10 mK bolometer using twenty 334g Te02 crystals to 
search for the decay [55]. 

Ionisation experiments : These passive experiments are mostly built in form 
of TPCs where the emitter is either the hlling gas or is included in thin foils. 
The advantage is that energy measurements and tracking of the two electrons 
is possible. Moreover, disadvantages are the energy resolution and the limited 
source strength by using thin foils. An experiment using a TPC with an active 
volume of 180 1 hlled with Xe (enriched to 62.5 % in which corresponds 

to 3.3 kg) under a pressure of 5 atm is done in the Gotthard-tunnel [56]. A 
TPG at UG Irvine was used to study ®^S'e , and ^^^Nd . A combina¬ 

tion of drift chambers, plastic scintillators and Nal-detectors is used in the 
ELEGANT V detector, investigating samples of the order of 100 g enriched 
in and ^^^Cd [57]. Enriched foils of , ^^^Cd and are 

also used by the NEMO-2 experiment [58]. 

Geochemical experiments : An alternative approach relies on the detection of 
the daughter nucleus. The geochemical method is using very old ores, which 
have accumulated a signihcant amount of daughter nuclei. Glearly the advan¬ 
tage of such experiments is the long exposure time of up to a billion years. 
However several new uncertainties are coming into consideration like an ac¬ 
curate age determination, to exclude other processes producing the daughter, 
avoid a high initial concentration of the daughter and to have a signihcant 
source strength. From all these considerations, only Se and Te-ores are usable. 
®^S'e , and decay to inert noble gases Xe) and the de¬ 

tection is based on isotopical anomalies due to double beta decay which are 
measured by mass spectrometry [59]. 

Radiochemical experiments : This method takes advantage of the radioactive 
decay of the daughter nuclei, allowing a shorter ’’measuring” time than geo¬ 
chemical experiments. They focus on the decay ‘^^‘^Th —jj 238^ _,.238 
Pu with Q-values of 850 keV and 1.15 MeV respectively. For the detection of 
the 238^ _i.238pg^ decay, the emission of a 5.5 MeV a-particle from the ^^®Pu 
decay is used as a signal. Of course geo- and radiochemical methods are not 
able to distinguish between the different double beta decay modes and are 
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finally limited in their sensitivity by 2h>(3(3 decay . 


2v(3(3 decay The predicted half-life for the 2v(3(3 decay is given by 

(T,%)-i = (39) 


where G'^'^ corresponds to the phase space and is the matrix element 

describing the transition. The main uncertainties in predicting accurate life¬ 
times are given by the errors on the matrix elements. A reliable knowledge 
of the matrix elements is necessary, because it influences the extractable neu¬ 
trino mass limit in the Qvj3j3 decay as well. In the past, it was quite common 
to work in the closure approximation, the replacement of the energies of the 
virtual intermediate states by an average energy, allowing the summation over 
all intermediate states because X) I 1= 1- Therefore only the wavefunc- 

tions of the initial and hnal state have to be known. But because interference 
between the different contributions has to be taken into account, all ampli¬ 
tudes have to be weighted with the correct energy and closure fails as a good 
description. The present determination of the matrix elements are done with 
QRPA-calculations. For details see [41,60-62]. All calculations depend on the 
strength of a particle-particle force which is a free parameter and has to be 
adjusted. A complete list of half-life calculations for A > 70 can be found in 
|63|. 

The hrst evidence for double beta decay was observed in geochemical experi¬ 
ments using selenium and tellurium-ores [64,65]. Newer measurements can be 
found in [59,66,67]. Because of phase space arguments, the detection of the 
decay has to be attributed to 2z//3/5 decay . A radiochemical detection 
of double beta decay using ^38^ with a half-life of 2.0 ±0.6 ■ 10^^ y [68] is con¬ 
sistent with theoretical predictions for 2vj3j3 decay . In 1987 the first direct 
laboratory detection by using ®^S'e was reported [69]. Meanwhile 2z//5/5 decay 
has been observed in several isotopes which are listed in Table 4. The highest 
statistics is obtained by the Heidelberg-Moscow experiment which has ob¬ 
served more than 20000 events (for comparison the first observation in 1987 
only consisted of 36 events). 


0z//5/5 decay The half-life for the Oz//?/? decay is given by (assuming rriu^ 
1 MeV) 



G^'^iMcT - Mpf 


me 


2 


(40) 
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Table 4 

Compilation of observed 2z//3/3 decay half-lives in several isotopes. * corresponds to 
geochemical measurements. 


Isotope 

Experiment 

Ti/2(1020y) 

48Ca 

Calt.-KIAE 

± 0.14 


MPIK-KIAE 

17.7 ±0.1+|;5 


IGEX 

11 ± 1.5 

82Se 

NEMO 2 

0.83 ±0.10 ±0.07 

ioomo 

ELEGANT V 

011 ^+*^-03 
U.iiO_o.o2 

ioomo 

NEMO 2 

0.095 ± 0.004 ± 0.009 

ioomo 

UCI 

0.0682lg;[j|j|± 0.0068 


NEMO 2 

0.375 ± 0.035 ± 0.021 


ELEGANT V 

'3-^0_o.o5 

128 T0* 

Wash. Uni-Tata 

77000 ± 4000 

150Nd 

ITEP/INR 

0.188;[|;|j^^ ± 0.019 

150Nd 

UCI 

0.0675;[j;[][j^^ ± 0.0068 


where the effective Majorana neutrino mass {niy^) is given by 

( 41 ) 

i 


with the relative CP-phases rji = ±1, Uei as the mixing matrix elements and 
mi as the corresponding mass eigenvalues. The expression can be generalised 
if right-handed currents are included to 


(TOh)-i = 

me 

+Cmri{- - ~){V) + - ~){^) + Cfjx{v)W 

me me 

where the coefficients C contain the phase space factors and the matrix ele¬ 
ments, 


iv) = vT.UejVe, {\) = XY^Ue.Ve, 

j 3 


(42) 


with Vej as the mixing matrix elements between right-handed neutrinos. 
Eq.(42) reduces to eq.(40) in case (77), (A) = 0. Also in 0z//3/5 decay the matrix 
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Table 5 

Compilation of neutrinoless double beta decay half-life and mass limits of the most 
investigated isotopes. The phase space factors and Q-values are taken from [47]. ^ 
after [73], * corresponds to geochemical measurement. 


Decay 

Q-value (keV) 

(G°")-'(y) 

nh (y) 

{rriy] 

> (eV) 

48(-( 

2ov^a 

-iTi 

4271 ± 4 

4.10E24 

> 9.5 • 1021(76%) 

< 

12.8(76%)^ 

76/-' 

32'-'6 


2039.6 ± 0.9 

4.09E25 

> 1.1 • 1025(90%) 

< 

0.5 

(90%) 

82c„_ 

34^16 

.iKr 

2995 ± 6 

9.27E24 

> 2.7-1022(68%) 

< 

5.0 

(68%) 



3034 ± 6 

5.70E24 

> 5.2 • 1022(68%) 

< 

5.0 

(68%) 

'ifcd- 

-'5'0'Sn 

2802 ± 4 

5.28E24 

> 2.9 • 1022(90%) 

< 

4.1 

(90%) 

'iTe- 

HifXe 

868 ± 4 

1.43E26 

> 7.7-1024(68%) 

< 

1.1 

(68%)* 

'f2°Te- 

HfOXe 

2533 ± 4 

5.89E24 

> 5.6 - 1022(90%) 

< 

3.0 

(90%) 

iflXe- 

136^ 

^ 56-oa 

2479 ± 8 

5.52E24 

> 4.4-1023(90%) 

< 

2.3 

(90%) 

1 

00 


3367.1 ± 2.2 

1.25E24 

> 2.1 - 1021(90%) 

< 

4.1 

(90%) 


element calculations are done with QRPA-calculations [63,70-72]. The gen¬ 
eral agreement between the calculations done by different groups are within a 
factor 2-3. 

From the experimental point, the evidence for 0z//5/5 decay in the sum energy 
spectrum of the electrons is a peak at the position corresponding to the Q- 
value of the involved transition. The half-life limits obtained so far for several 
different isotopes are shown in Table 5. The best limit is coming from the 
Heidelberg-Moscow experiment resulting in a bound of [74] (Fig.5) 

T0/2 > 1-1 ■ ^ < 0.47614 (90%C'L) (43) 

using the matrix elements of [63]. Because in most see-saw models 
corresponds to [10] (see Table 1), this bound is much stronger than 
single beta decay but applies only if neutrinos are Majorana particles. Al¬ 
lowing also right-handed currents to contribute, is hxed by an ellipsoid 

which is shown in Fig. 6. As can be seen, the largest mass allowed occurs for 
(A), (77) 7^ 0. In this case the half-life of eq. 43 corresponds to 


< 0.56e14 

(44) 

{ri) < 6.5 -10"® 

(45) 

(A) < 8.2 -10"^ 

(46) 


The limit also allows a bound on a possible right-handed Wr which is shown 
in Fig. 7. Together with vacuum stability arguments a mass for Wr lower 
than about 1 TeV can be excluded. The influence of double charged Higgses, 
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Fig. 5. Observed spectrum near the expected OulSP decay peak of the Heidel- 
berg-Moscow collaboration. No signal can be seen. The two different spectra corre¬ 
spond to measuring periods with (black) and without (grey) pulse shape discrimi¬ 
nation (with kind permission of H.V. Klapdor-Kleingrothaus). 



Fig. 6. Obtained neutrino mass {my^) as a function of the two right-handed current 
parameters (A), {rj). On the right the projection on the (A), {rf) plane is shown. 

which also can contribute to neutrinoless double beta decay is shown as well 
[75]. Limits on other interesting quantities like the R-parity violating SUSY 
parameter [76] and leptoquarks [77] can be derived. From the point of 
right-handed currents the investigation of the transition to the first excited 
state is important, because the mass term here vanishes in first order. The 
phase space for this transition is smaller, but the de-excitation photon might 
allow a good experimental signal. For a compilation of existing bounds on 
transitions to excited states see [78]. As long as no signal is seen, bounds from 
ground state transitions are much more stringent on right-handed parameters. 
Eq.(41) has to be modified in case of heavy neutrinos {rriy^ 1 MeV). For 
such heavy neutrinos the mass can no longer be neglected in the neutrino 
propagator resulting in an A-dependent contribution 

N M 

=1 H H F{mh,A)U^f^mh\ (47) 

i=l,light h=l,heavy 
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Fig. 7. Bound on the mass of a right-handed W as a function of the heavy neutrino 
mass ruAT. under the assumption of different masses of the double charged higgs 
A . Shown are the regions excluded by double beta decay (lower left) and from 
vacuum stability (lower right). The curves a-e correspond to masses of 0.3 TeV, 1 
TeV, 2 TeV, 5 TeV and infinite for the A (from [75]). 


By comparing these limits for isotopes with different atomic mass, interesting 
limits on the mixing angles and v-r parameters for an MeV v-r can be obtained 
[79,80]. 

A complete new class of decays emerges in connection with majoron emis¬ 
sion in double beta decay [81]. The majoron y is the Goldstone-boson of a 
spontaneous breaking of a global lepton-number symmetry. Depending on its 
transformation properties under weak isospin, singlet [82], doublet [83] and 
triplet [84] models exist. The triplet and pure doublet model are excluded 
by the measurements of the Z-width at LEP, because they would contribute 
2 (triplet) or 0.5 (doublet) neutrino flavours. Several new majoron-models 
evolved during the last years [85,86]. In consequence different spectral shapes 
for the sum electron spectrum are predicted which can be written as 

dN 

^{Q-Er-F{E,Z) (48) 


where F(E,Z) is the Fermi-function and the spectral index n is 1 for the clas¬ 
sical majoron, n=3 for lepton number carrying majorons, n=5 for decay 
and n=7 for several other majoron models. A different shape is obtained in the 
vector majoron picture of Carone [87]. It should be noted that supersymmetric 
Zino-exchange allows the emission of two majorons, which also results in n=3, 
but a possible bound on a Zino-mass is less stringent than direct accelerator 
experiments [88]. In the n=l model the effective neutrino -majoron coupling 
(Qux) can be deduced from 

(T°/7)-i =1 Mgt - Mf |2 I |2 (49) 
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Table 6 

Compilation of limits on half-lives and effective majoron-neutrino couplings from 
different isotopes. * corresponds to a geochemical measurement, which has no power 
to discriminate between different decay modes. 


Isotope 

Experiment 

Ai/2(102iy) 

(5.x)(10-") 

48Ca 

ITEP 

0.72 (90 %) 

5.3 


MPIK-KIAE 

7.91 (90 %) 

2.3 


ITEP 

10 (68 %) 

2.2 

^^^Ge 

UCSB-LBL 

1.4 (90 %) 

5.8 

76Ge 

PNL-USC 

6.0 

2.8 


Cal.-PSI-Neu 

1.0 (90 %) 

6.9 

82Se 

NEMO 2 

2.4 (90 %) 

1.4 

ioomo 

ELEGANT V 

5.4 (68 %) 

0.7 

ioomo 

NEMO 2 

0.5 (90 %) 

2.3 

ioomo 

UCI 

0.3 (90 %) 

3 

128 T0* 

Wash. Uni-Tata 

7700 

0.3 

i36Xe 

Cal.-PSI-Neu 

14 (90 %) 

1.5 

150Nd 

UCI 

0.28 (90 %) 

1 


where is given by 

= 'T.duxUeiUej (50) 

hi 


Present half-life limits for this decay (n=l) and the deduced coupling constants 
are given in Table 6. A hrst half-life limit for the n=3 mode was given in 
[89], a evaluation for is given in [90]. A more recent extensive study 

of all modes can be found in [86]. Limits obtained by the Heidelberg-Moscow 
experiment with a statistical signihcance of 4.84 kg-y are [91] 


> 5.85 ■ lO^^i/ (n=3) (51) 

> 6.64 ■ lO^^y (n=7) (52) 

Additionally the j3^j3^-(\.ec&y in combination with EC can be observed via 
the following decay modes 


(Z, A)^{Z- 2, A) + 2e+ + (2z/e) (53) 

A i^i ^ ~ 2, A) -|- e"*" -|- (2z/e) (54) 
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(55) 

( 56 ) 


2ej5 + (^Z, A) —> (^Z — 2, A) + 2z/e 
2e^ + {Z, A) ^ {Z — 2, A)* —{Z — 2, A) + 'j + 2X-rays 

is always accompanied by EC/EC or /5+/EC-decay. Because of the 
Coulomb-barrier and the reduction of the Q-value by 4 rrieC^, the rate for 
is small and energetically only possible for seven nuclides. Predicted half-lives 
for j3^j3^ are of the order 10^® y while for /d^/EC this can be reduced by orders 
of magnitude down to y making an experimental detection more real¬ 

istic. The experimental signature of the decay modes is rather clear because of 
the two or four 511 keV photons. The last mode (eq.56) to an excited state is 
giving a characteristic gamma associated with X-ray emission. Half-lives ob¬ 
tained with ^'^®Cd and '^^Ru are of the order 10^® y [92,78]. Extracted neutrino 
mass limits are orders of magnitude worse than the Qvj3j3 decay limits, but if 
there is any positive observation of the Oz//?/? decay mode, the /5’''/EC-mode 
can be used to distinguish whether this is dominated by the neutrino mass 
mechanism or right-handed currents [93]. 


Future Several upgrades are planned to improve the existing half-life limits. 
Because of the enormous source strength after additional years of running 
the dominant project will still be the Heidelberg-Moscow experiment probing 
neutrino masses down to 0.2 eV. A new experiment to improve the sensitivity 
on is ELEGANT VI, using 25 modules of CaF 2 with a total amount of 
31 g of within a Csl detector [57]. A different approach might be the 
use of CaF 2 {Eu) as a cryogenic bolometer and to measure simultaneously 
the scintillation light [94]. is interesting because it can be treated with 
nuclear shell model calculations. The building up of NEMO-3, which should 
start operation in 1999, will allow to use up to 10 kg of material in form 
of foils for several isotopes like [95]. Even more ambitious would be 

the usage of large amounts of materials (in the order of several hundred kg 
to tons) like enriched added to scintillators [96], 750 kg Te02 in form 

of cryogenic bolometers (CHORE) [97] or a huge cryostat containing several 
hundred detectors of enriched with a total mass of 1 ton (GENIUS) [98]. 
Further, ideas to use a large amount of and detect the created daughter 

^^®Ba with atomic traps and resonance ionisation spectroscopy exist. This will 
allow no information on the decay mode and will be dominated by decay 
[99,100]. 

4-2 Magnetic moment of the neutrino 


Another possibility to check the neutrino character is the search for its mag¬ 
netic moment . In the present standard model both types of neutrinos have 
no magnetic moment because neutrinos are massless and a magnetic moment 
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would require a coupling of a left-handed state with a right-handed one which 
is absent. A simple extension by including right-handed singlets allows for 
Dirac-masses. In this case, it can be shown that neutrinos can get a magnetic 
moment due to loop diagrams which is proportional to their mass and is given 
by [101,102] 

“ sTS?’""" ' * 

In case of neutrino masses in the eV-range, this is far to small to be observed 
and to have any signihcant effects in astrophysics. Nevertheless there exist 
GUT-models, which are able to increase the magnetic moment without in¬ 
creasing the mass [103]. However Majorana neutrinos still have a vanishing 
static moment because of CPT-invariance. The existence of diagonal terms 
in the magnetic moment matrix would therefore prove the Dirac-character of 
neutrinos . Non-diagonal terms in the moment matrix are possible for both 
types of neutrinos allowing transition moments of the form z/g - 
Limits on magnetic moments arise from z/e e - scattering experiments and as- 
trophysical considerations. The differential cross section for e - scattering 
in presence of a magnetic moment is given by 


^ + X + QaY + {Qv + X - 

2.meT^ ^ Tia'^IJ,Y ^ - T/E„ 


+{9l-{x + 9vy)-Yr] + 


Ey ' ml T 
where T is the kinetic energy of the recoiling electron and 


gv = 2siy9w + - 9a = ±- (+(-) for yy) 


(58) 

(59) 


(60) 


and X denotes the neutrino form factor related to its square charge radius (r^) 
X = — -^y)siy9w X — >• —X for z7g (61) 

The contribution associated with the charge radius can be neglected in the 
case fi„ ~ As can be seen, the largest effect of a magnetic moment 

can be observed in the low energy region, and because of destructive interfer¬ 
ence of the electroweak terms, searches with antineutrinos would be preferred. 
The obvious sources are therefore nuclear reactors. Experiments done so far 
result in a bound of < 1.52 ■ 10“^Vs for z/g [104]. Measurements based on 
z/g e —z/gC and z/^ e —> z/^e scattering were done at LAMPF and BNL yielding 
bounds for z/g and z/^ of ~ 10“®/iB [105,106]. 
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Astrophysical limits are somewhat more stringent but also more model de¬ 
pendent. An explanation of the solar neutrino problem by spin precession of 
vl into vn done by the magnetic field of the solar convection zone requires 
a magnetic moment of the order k, [107]. Observation of 

neutrinos from Supernova 1987A yield a somewhat model dependent bound of 
Hu < [108,109]. Also the neutrino emissivity of globular cluster stars 

done by excessive plasmon decay 7 —z/z/ is only consistent with observation 
for a magnetic moment of the same order [110]. This last bound applies to 
neutrinos lighter than 5 keV. 

To improve the experimental situation and especially check the region relevant 
for the solar neutrino problem new experiments are under construction. The 
most advanced is the NUMU experiment [111] currently installed near the 
Bugey reactor. It consists of a 1 m^ TPC loaded with CF 4 under a pressure 
of 5 bar. The usage of a TPC will not only allow to measure the electron 
energy but for the first time in such experiments also the scattering angle, 
therefore allowing the reconstruction of the neutrino energy. The neutrino en¬ 
ergy spectrum at reactors in the energy region 1.5 < 77,^ < 8 MeV is known 
at the 3 % level. To suppress background, the TPC is surrounded by 50 cm 
anti-Compton scintillation detectors as well as a passive shielding of lead and 
polyethylene. In case of no magnetic moment the expected count rate is 9.5 per 
day increasing to 13.4 per day if Hv = for an energy threshold of 500 

keV. The estimated background is 6 events per day. The expected sensitivity 
level is down to Hu = ^ ■ 10“^^/iB . The usage of a low background Ge-Nal 
spectrometer in a shallow depth near a reactor has also been considered [112]. 
The usage of large low-level detectors with a low-energy threshold of a few keV 
in underground laboratories is also under investigation. The reactor would be 
replaced by a strong /^-source. Calculations for a scenario of a 1-5 MCi ^"^^Pm 
source (endpoint energy of 234.7 keV) in combination with a 100 kg low-level 
Nal(Tl) detector with a threshold of about 2 keV can be found in [113]. 


4-3 Search for heavy Majorana neutrinos 


For the see-saw-mechanism to work heavy Majorana neutrinos N are nec¬ 
essary. The required lightness of neutrino masses makes a detection of the 
corresponding heavy state impossible. The mixing of a heavy neutrino mu to 
the light state is ruled by 6* = ^JmL/mH 1- However there exist the¬ 
oretical models which decouple the mixing from any mass relation [114,115]. 
Assuming that in eq.(3) 7^ 0 and that by an internal symmetry at tree 

level the relation mffmff = is valid, the mixing is decoupled from the 

ratio mi/m 2 and can be close to one in case that mff ~ bt//. Masses for light 
neutrinos vanish at tree level and will be generated at higher orders. 

From the experimental point of view, heavy Majorana neutrinos can be 
searched for at accelerators. The LEP-data on the Z°-width already exclude 
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any additional neutrino lighter than 45 GeV. Searches for heavier neutrinos 
have been done at LEP1.5. The search for Majorana neutrinos heavier than 
the Z° focusses on the N-decay channels 

and N ^ uZ° (62) 

which is identical to signatures looked for in searches of excited fermions. A 
detailed description of pair production of heavy Dirac and Majorana neutrinos 
in e’''e“collisions can be found in [116]. Pair production of Majorana neutrinos 
would result in two like-sign charged leptons. Furthermore, HERA offers the 
chance to search for heavy Majorana masses in ep-collisions [117]. For accu¬ 
mulated 200 ph~^ a discovery limit up to 160 GeV is possible. Also future 
high energy machines allow an extended search for heavy neutrinos via 
reactions 


e e 


vN N ^ ^ vZ,N ^ vH 


(63) 


The dominant background will be lV’''fP“production [118]. LHG offers 
searches either in the pair-production or single Majorana neutrino produc¬ 
tion mode [119-121]. The advantage of single Majorana production is that it 
depends only linearly on the neutrino mixing. The single production channel 
via 


pp e-NX e-e-fP+X 


(64) 


offers a signal of two same sign leptons, two jets with the invariant mass 
of and no missing energy. For an assumed luminosity of 10 fb“^ the 
discovery potential goes up to 1.4 TeV (0.8 TeV) for an assumed mixing of 
sin9 ^ 10“^(10“^). 


5 Neutrino oscillations 


In case of massive neutrinos the mass eigenstates do not have to be identical 
with the flavour eigenstates, similar to the GKM-mixing in the quark sector. 
This offers the chance for neutrino oscillations. Oscillations might be the only 
chance to see effects of and Ur in the eV mass range which is not accessible 
in direct experiments. 
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5.1 General formalism 


The concept of neutrino oscillations has been introduced by [122], The weak 
eigenstates are related to the mass eigenstates z/j via a unitary matrix U 

H UaiVi (65) 

i 


which is given for Dirac neutrinos as 



C12C13 

-S12C13 

Si3e 


iS 

iS 


—"S12C23 

- Ci2S23Sl3e“ 

C12C23 - -Si2S23-Sl3e“ 

"S23C13 


iS 

iS 


'Sl2'S23 ■ 

- Ci2S23Sl3e“ 

-C12S23 - Si2C23Sl3e 

C23C13 


( 66 ) 


and in the Major ana case as 

( 

U = 

y 5125236"^''"" ' - Ci2S23Si3e 


C 12 C 13 

Sl2C23e*^^" - Ci2S23Sl3e*^^i®+^"3) 


-Sl2Cl3e 

C12C23 - Si2S23Sl3e*^^"®+^i®"'^i"^ 

-Cl2S23e'^^^ - Si2C23Si3e'('5«-5i2) 


Si3e-*'^13 ^ 

S23Cl3e*‘^23 


(67) 


C 23 C 13 j 


with c,s = cos6, sinO. The quantum mechanical transition probability can be 
derived (assuming relativistic neutrinos and CP-conservation) as [7] 

P(z/„ ^ z/^) = ^ I U/si n Uo,i P +ReJ2 U0iU0jU*iUajexp{-itAm‘^ij/2E){68) 
i 


with Am^ij =1 — m'j \. In the simple two-flavour mixing the probability to 

hnd Up in a distance x with respect to a source of z/q, is given by 

TTX 

P(z/q, — Up) = sin^29sin^ — (69) 

Jj 


giving the oscillation length L in practical units as 


AnEh 

Am^G' 


2.48( 


E 

MeV 


eV 

Am^ 


m 


(70) 


For a more extensive review on N flavour mixing, wave-packet treatment and 
coherence considerations see [7,123,124], Terrestrial experiments are done with 
nuclear reactors and accelerators. The discussed oscillation searches involve 
the three known neutrinos as well as a possible sterile neutrino z/ 5 . 
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5.2 Reactor experiments 


Reactor experiments are disappearance experiments looking for Vf. ^ vx- 


5.2.1 Principles 

Reactors are a source of MeV z/g due to the fission of nuclear fuel. The main 
isotopes involved are and ^^^Pu. The neutrino rate per hssion 

has been measured [126] for all isotopes except ^sspj jg good agreement 
with theoretical calculations [127]. Experiments typically try to measure the 
positron spectrum which can be deduced from the z/g - spectrum and either 
compare it directly to the theoretical predictions or measure it at several 
distances from the reactor and search for spectral changes. Both types of 
experiments were done in the past. The z/g cross section is known to about 1.4 
% [128]. The detection reaction is 

^e + p^e'^ + n (71) 


with an energy threshold of 1.804 MeV. The detection reaction (71) is always 
the same, resulting in different strategies for the detection of the positron and 
the neutron. Normally coincidence techniques are used between the annihi¬ 
lation photons and the neutrons which diffuse and thermalise within 10-100 
/iS. The main background are cosmic ray muons producing neutrons in the 
surrounding of the detector. 


5.2.2 Experimental status 

Several reactor experiments have been done in the past (see Table 7). All 
these experiments had a hducial mass of less than 0.5 t and the distance to 
the reactor was never more than 250 m. Two new reactor experiments pro¬ 
ducing data are CHOOZ and Palo Verde. The CHOOZ-experiment in France 
[129] has some advantages with respect to previous experiments. First of all 
the detector is located underground with a shielding of 300 mwe, reducing 
the background due to atmospheric muons by a factor of 300. Moreover, the 
detector is about 1030 m away from two 4.2 GW reactors (more than a factor 
4 in comparison to previous experiments) enlarging the sensitivity to smaller 
Am^ . In addition the main target has about 4.8 t and is therefore much larger 
than those used before. The main target consists of a specially developed Gd- 
loaded scintillator. This inner detector is surrounded by an additional detector 
containing 17 t of scintillator without Gd and 90 t of scintillator as an outer 
veto. The signal in the inner detector is the detection of the annihilation pho¬ 
tons in coincidence with n-capture on Gd, the latter producing gammas with 
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Table 7 

List of finished reactor experiments. Given are the power of the reactors and the 
distance of the experiments with respect to the reactor. 


reactor 

thermal power [MW] 

distance [m] 

ILL-Grenoble (F) 

57 

8.75 

Bugey (F) 

2800 

13.6, 18.3 

Rovno (USSR) 

1400 

18.0,25.0 

Savannah River (USA) 

2300 

18.5,23.8 

Gosgen (CH) 

2800 

37.9, 45.9, 64.7 

Krasnojarsk (Russia) 

? 

57.0, 57.6, 231.4 

Bugey III (F) 

2800 

15.0, 40.0, 95.0 
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Fig. 8. The measured positron spectrum of the CHOOZ-experiment (left). On the 
right the ratio with the expectation is shown. No oscillation signal is visible (from 
[130]). 

a total sum of up to 8 MeV. The first published positron spectrum [130] is 
shown in Fig. 8 and shows no hints for oscillations. The resulting exclusion 
plot is shown in Fig.18. 

The second experiment is the Palo Verde (former San Onofre) experiment [131] 
near Phoenix, AZ (USA). It consists of 12 t liquid scintillator also loaded with 
Gd. The scintillator is hlled in 66 modules arranged in an 11x6 array. The 
coincidence of three modules serves as a signal. The experiment is located un¬ 
der a shielding of 46 mwe in a distance of about 750 (820) m to three reactors 
with a total power of 10.2 GW. 

A further project plans a 1000 t liquid scintillator detector (KamLAND) [132]. 
It is approved by the Japanese Government and will be constructed at the 
Kamioka site. Having a distance of 160 km to the next reactor, it will probe 
Am^ down to 10“^eU^. 
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Fig. 9. Neutrino spectrum from 7r-decay. Beside the monoenergetic line from 
7r+-decay at rest, the continuous spectra of v^. and are shown. 

5.3 Accelerator experiments 


The second source of terrestrial neutrinos are high energy accelerators. Exper¬ 
iments can be of either appearance or disappearance type [125]. 


5.3.1 Principles 

Accelerators typically produce neutrino beams by shooting a proton beam on 
a hxed target. The produced secondary pions and kaons decay and create a 
neutrino beam dominantly consisting of . The detection mechanism is via 
charged weak currents 


i>iN i + X i = e, 


(72) 


where N is a nucleon and X the hadronic hnal state. Depending on the intended 
goal, the search for oscillations therefore requires a detector which is capable 
of detecting electrons, muons and r - leptons in the hnal state. 


5.3.2 Experimental status 

Accelerators at medium energy At present there are two experiments 
running with neutrinos at medium energies {Ei, 30 - 50 MeV) namely KAR¬ 
MEN and LSND . Both experiments use 800 MeV proton beams on a beam 
dump to produce pions. The expected neutrino spectrum from pion and p- 
decay is shown in Fig. 9. The beam contamination of Ug is in the order of 
10“^. The KARMEN experiment [133] at the neutron spallation source ISIS 
at Rutherford Appleton Laboratory is using 56 t of a segmented liquid scintil¬ 
lator. The main advantage of this experiment is the known time structure of 
the two proton pulses hitting the beam dump (two pulses of 100 ns with a sep- 
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Fig. 10. Vf. - parameter plot. Shown is the region of evidence from the LSND 
vr-decay at rest analysis (grey areas) and different exclusion curves from Bugey, 
KARMEN, E776 (left) and the preliminary NOMAD result (right). 

aration of 330 ns and a repetition rate of 50 Hz). Because of the pulsed beam, 
positrons are expected within 0.5-10.5 ps after beam on target. The signature 
for detection is a delayed coincidence of a positron in the 10 - 50 MeV region 
together with 7 -emission from either p(n, 7 )D or Gd(n, 7 )Gd reactions. The 
hrst results in 2.2 MeV photons while the latter allows gammas up to 8 MeV. 
The limit reached so far is shown in Fig. 10. Recently KARMEN published 
a 2- and 3- flavour analysis of z/g -Vt and z/g - z/^ oscillations by comparing 
the energy averaged GG-cross section for interactions with expectation as 
well as making a detailed maximum likelihood analysis of the spectral shape 
of the electron spectrum observed from (z/g ,e“)^^V gg reactions [134]. To 
improve the sensitivity by reducing the neutron background, a new veto shield 
against atmospheric muons was constructed which has been in operation since 
Feb. 1997 and is surrounding the whole detector. The region which can be in¬ 
vestigated in 2-3 years of running in the upgraded version is also shown in 
Fig. 10. 

The LSND experiment [135] at LAMPF is a 167 t mineral oil based liquid 
scintillation detector using scintillation and Gerenkov light for detection. It 
consists of an approximately cylindrical tank 8.3 m long and 5.7 m in diame¬ 
ter. The experiment is about 30 m away from a copper beam stop under an 
angle of 12° with respect to the proton beam. For the oscillation search in 
the channel Vg -Pg a signature of a positron within the energy range 36 MeV 
< Ee < 60 MeV together with an in time and spatial correlated 2.2 MeV pho¬ 
ton from p(n, 7 )D is required. The analysis (Fig. 11) [136] ends up in evidence 
for oscillations in the region shown in Fig. 10. Recently LSND published their 
Vf. - Vg analysis for pion decays in flight by looking for isolated electrons in the 
region 60 MeV < < 200 MeV coming from (z/g ,e“)^^V gg reactions 

[137], which is in agreement with the former evidence from pion decay at rest. 
Also LSND continues with data acquisition. 

An increase in sensitivity in the Vg - v^. oscillation channel can be reached in 
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Fig. 11. Positron spectrum of beam-excess data. The dashed curve corresponds to 
the estimated neutrino background while the solid line indicates neutrino oscillations 
for large Am? plus the estimated background (from [136]). 

the future if there is a possibility for neutrino physics at the planned European 
Spallation Source (ESS) or the National Spallation Neutron Source (NSNS) at 
Oak Ridge which might have a 1 GeV proton beam in 2004. The Fermilab 8 
GeV proton booster offers the chance for a neutrino experiment as well which 
could start data taking in 2001. It would use part of the LSND equipment and 
consist of 600 t mineral oil contained and be located 500 m away from the 
neutrino source (MiniBooNE)[138]. An extension using a second detector at 
1000m is possible (BooNE). At GERN the PS neutrino beam could be revived 
with an average energy of 1.5 GeV and two detector locations at 128 m and 850 
m as it was used by the former GDHS [139] and GHARM-experiment [140]. 
By measuring the Uf, / ratio the complete LSND region can be investigated 
[141], 


Accelerators at high energy High energy accelerators provide neutrino 
beams with an average energy in the GeV region. With respect to high en¬ 
ergy experiments at present especially GHORUS and NOMAD at GERN will 
provide new limits. They are running at the GERN wide band neutrino beam 
with an average energy of around 25 GeV, produced by 450 GeV protons ac¬ 
celerated in the SPS and then hitting a beryllium beam dump. To reduce the 
uncertainties in the neutrino flux predictions, the NA56 - experiment mea¬ 
sured the resulting pion and kaon spectra [142]. The experiments are 823 m 
(GHORUS) and 835 m (NOMAD) away from the beam dump and designed 
to improve the existing limits on - Uj- oscillations by an order of magnitude. 
The beam contamination of prompt from D^-decays is of the order 2-5 
■10“® [143,144]. Both experiments differ in their detection technique. While 
GHORUS relies on seeing the track of the r - lepton and the associated decay 
vertex with the kink because of the r-decay, NOMAD relies on kinematical 
criteria. 

The GHORUS experiment [145] uses emulsions with a total mass of 800 kg 
segmented into 4 stacks, 8 sectors each as a main target. To determine the 
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vertex within the emulsion as accurate as possible, systems of thin emulsion 
sheets and scintillating fibre trackers are used. Behind the tracking devices 
follows a hexagonal air core magnet for momentum determination of hadronic 
tracks, an electromagnetic lead-scintillating fibre calorimeter with an energy 
resolution of AE/E = 13%/\/E for electrons as well as a muon spectrometer. 
At- lepton created in the emulsion by a charged current reaction is pro¬ 
ducing a track of a few hundred pm. After the running period the emulsions 
are scanned with automatic microscopes coupled to CCDs. The experiment 
searches for the muonic and hadronic decay modes of the r and took data 
from 1994 to 1997. The present limit (Fig. 12) provided by CHORUS for the 
Ufj_ - Ur channel for large Am^ is [146] 

sin‘^2e < 1.2 X 10"^ (90%C'L) (73) 


The final goal is to reach a sensitivity down to siv?29 k, 2 x 10“^ for large 
Arri^ . 

The NOMAD experiment [147] on the other hand relies on the kinematics. 
It has as a main active target 45 drift chambers representing a total mass 
of 2.7 tons followed by transition radiation and preshower detectors for e/vr 
separation. After an electromagnetic calorimeter with an energy resolution of 
AE/E = 3.22 %/^/E © 1.04% and a hadronic calorimeter five muon chambers 
follow. Because most of the devices are located within a magnetic field of 0.4 
T a precise momentum determination due to the curvature of tracks is possi¬ 
ble. The r-lepton cannot be seen directly, the signature is determined by the 
decay kinematics. The main background for the r-search are regular charged 
and neutral current reactions. In normal charged current events the muon 
balances the hadronic final state in transverse momentum with respect 
to the neutrino beam. Hence the value for missing transverse momentum is 
small. The angle ^ih between the outgoing lepton and the hadronic final state 
is close to 180" while the angle ^mh between the missing momentum and the 
hadronic final state is more or less equally distributed. In case of a r - decay 
there is significant missing px because of the escaping neutrinos as well as a 
concentration of ^mh to larger angles because of the kinematics. In the u^ - 
Ur channel for large Am? NOMAD gives a limit of [148] 

si7?2e < 1.2 X 10-3 (90%C'L) (74) 


Both limits are now better than the limit of E531 (Fig.12). Having a good 
electron identification, NOMAD also offers the possibility to search for os¬ 
cillations in the u^ - Ue channel. A preliminary limit (Fig.10) on - z/g is 
available as (for large Am?) 

sir?2e < 2 X 10-3 (90%C'L) (75) 
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This and a recently published CCFR result [149] seem to rule out the large 
Am^ region of the LSND evidence. NOMAD will continue data taking until 
Sept. 1998. 


5.3.3 Future experiments 

Possible future ideas split into two groups depending on the physics goal. 
One group is focussing on improving the existing bounds in the eV-region 
by another order of magnitude with respect to CHORUS and NOMAD. This 
effort is motivated by the classical see-saw-mechanism which offers a Ur in the 
eV-region as a good candidate for hot dark matter by assuming that the solar 
neutrino problem can be solved by z/g - oscillations. The second motivation 
is to check the LSND evidence. The other group plans to increase the source 
- detector distance to probe smaller Am^ and to be directly comparable to 
atmospheric scales (see chapter 5.4). 


Short and medium baseline experiments Several ideas exist for a next 
generation of short and medium baseline experiments. At CERN the proposed 
follow up is TOSCA, a detector combining features of NOMAD and CHORUS 
[150]. The idea is to use 2.4 tons of emulsion within the NOMAD magnet in 
form of 6 target modules. Each module contains an emulsion target consisting 
of 72 emulsion plates, as well as a set of large silicon microstrip detector planes 
and a set of honeycomb tracker planes. Both will allow a precise determination 
of the interaction vertex improving signihcantly the efficiency. To verify the 
feasibility of large silicon detector planes maintaining excellent spatial reso¬ 
lution over larger areas, NOMAD included a prototype (STAR) in the 1997 
data taking. Moreover options to extract a neutrino beam at lower energy 
of the proton beam (350 GeV) at the CERN SPS to reduce the prompt Ur 
background are discussed. The proposed sensitivity in the - z/,- channel is 
around 2x10“^ for large Am^ (Am^ > lOOeU^) (Fig. 12). Also proposals for 
a medium baseline search exist [151,152]. The CERN neutrino beam used by 
CHORUS and NOMAD is coming up to the surface again in a distance of 
about 17 km away from the beam dump. An installation of an ICARUS-type 
detector (liquid Ar TPC) [151] could be made here. In a smaller version, two 
hne grained calorimeters located at CERN and in 17 km distance might be 
used as well [152]. 


Long baseline experiments Several accelerators and underground lab¬ 
oratories around the world offer the possibility to perform long baseline ex¬ 
periments . This is of special importance to probe the region of atmospheric 
neutrinos directly. 
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(present) CHORUS ^ NOMAD (present) 



Fig. 12. Status and proposed curves for - Ur oscillations. Shaded are the regions 
expected for an eV-Ur as dark matter candidate according to the see-saw-model, 
where Ue - oscillations are used to explain the solar neutrino deficit. Meanwhile 
the CHORUS and NOMAD exclusion curve moved down to 1.2 -10“^ for large Am? 


KEK - Super-Kamiokande : The first of these experiments will be the KEK- 
E362 experiment (K2K) [153] in Japan sending a nentrino beam from KEK 
to Snper-Kamiokande . The distance is 235 km. A 1 kt near detector, abont 
1 km away from the beam dnmp will serve as a reference and measnre the 
nentrino spectrum. The neutrino beam with an average energy of 1 GeV is 
produced by a 12 GeV proton beam dump. The detection method within 
Snper-Kamiokande will be identical to that of their atmospheric neutrino de¬ 
tection. The beamline should be hnished by the end of 1998 and the experiment 
will start data taking in 1999. In connection with the JHG-project an upgrade 
of KEK is planned to a 50 GeV proton beam, which could start producing 
data around 2004 and would make a Vr appearance search possible. 

Fermilab - Soudan: A neutrino program is also associated with the new Main 
Injector at Fermilab. The long baseline project will send a neutrino beam to 
the Soudan mine about 735 km away from Fermilab. Here the MINOS ex¬ 
periment [154] will be installed. It also consists of a near detector located at 
Fermilab and a far detector at Soudan. The far detector will be made of 8 
kt magnetized Fe toroids in 600 layers with 2.54 cm thickness interrupted by 
about 32000 m^ active detector planes in form of plastic scintillator strips 
with X and y readout to get the necessary tracking informations. An addi¬ 
tional hybrid emulsion detector for r-appearance is also under consideration. 
The project could start at the beginning of next century. 

CERN - Gran Sasso: A further program considered in Europe are long base¬ 
line experiments using a neutrino beam from GERN down to Gran Sasso 
Laboratory. The distance is 732 km. Several experiments have been proposed 
for the oscillation search. The hrst proposal is the IGARUS experiment [155] 
which will be installed in Gran Sasso anyway for the search of proton decay 
and solar neutrinos. This liquid Ar TPG can also be used for long baseline 
searches. A prototype of 600 t is approved for installation which will happen 
in 1999. A second proposal, the NOE experiment [156], plans to build a giant 
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lead-scintillating fibre detector with a total mass of 4.3 kt. The calorimeter 
modules will be interleaved with transition radiation detectors with a total of 
2.4 kt. The complete detector will have twelve modules, each 8mx8mx5m, and 
a module for muon identihcation at the end. A third proposal is the building of 
a 125 kt water-RICH detector (AQUA-RICH) [157], which could be installed 
outside the Gran Sasso tunnel. The readout will be done by 3600 HPDs with 
a diameter of 250 mm and having single photon sensitivity. Finally there ex¬ 
ists a proposal for a 750t iron-emulsion sandwich detector (OPERA) [158] 
which could be installed either at the Fermilab-Soudan or the CERN-Gran 
Sasso project. It could consist of 92 modules, each would have a dimension 
orthogonal to the beam of 3x3 m^ and would consist out of 30 sandwiches. 
One sandwich is composed out of 1 mm iron, followed by two 50 /rm emulsion 
sheets, spaced by 100 /rm. After a gap of 2.5 mm, which could be hlled by 
low density material, two additional emulsion sheets are installed. The r, pro¬ 
duced by CC reaction in the iron, decays in the gap region, and the emulsion 
sheets are used to verify the kink of the decay. 

A project in the very far future could be oscillation experiments involving 
a /i'''/i“-collider currently under investigation. The created neutrino beam is 
basically free of Vr and can be precisely determined to be 50 % and 

50% z/e iye) for Because the /r+/r“-collider would be a high luminos¬ 

ity machine, one even can envisage very long baseline experiments e.g. from 
Fermilab to Gran Sasso with a distance of 9900 km [159]. 


5-4 Atmospheric neutrinos 


A different source of neutrinos are cosmic ray interactions within the atmos¬ 
phere. A detailed prediction of the expected flux depends on three main in¬ 
gredients, namely the cosmic ray spectrum and composition, the geomagnetic 
cutoff and the neutrino yield of the hadronic interaction in the atmosphere. At 
lower energies (i7~ 1 GeV) neutrinos basically result from pion- and muon- 
decay leading to rough expectations for the fluxes like r'/i ~ ~ 2z/e or 

z/e/t'e ~ The ratio z/g / z/^ drops quickly above 1 GeV, because for 

Efj_ > 2.5 GeV the path length for muons becomes larger than the produc¬ 
tion height. At even higher energies the main source of z/g are A'es-decays 
(Kl —TTCz/e). Gontributions of prompt neutrinos from charm decay are neg¬ 
ligible and might become important in the region above 10 TeV [160]. For 
neutrinos in the energy range 300 MeV E, ^ 3 GeV the energy of the pri¬ 
mary typically lies in the region 5 GeV ^ Ep ~ 50 GeV. This region of the 
spectrum is affected by the geomagnetic cutoff, which depends on the gyro- 
radius of the particles, introducing a factor A/Z between nuclei and protons 
of the same energy. However neutrino production depends on the energy per 
nucleon E/N. Furthermore the energy range below about 20 GeV is also af¬ 
fected by the 11-year activity cycle of the sun, which is in the maximum phase 
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Muon — neutr'ino Response 



Fig. 13. Different classes of observable muons as a function of neutrino energy. The 
peak corresponding to contained events is reduced by a factor 10 (after [262]). 

preventing low energy cosmic rays to penetrate into the inner solar system. 
Neutrinos in the region well beyond 1 GeV can be detected via horizontal or 
upward going muons produced by CC reactions. The dominant part is given 
from events between 10 GeV < Ey < 10*^ GeV. The contribution of primaries 
with energies larger than 10^ GeV/nucleon to the upward going muon flux 
is only about 15 %. Several authors made calculations for the neutrino flux 
for different detector sites covering the energy region from 100 MeV to 10^ 
GeV [161-166]. The absolute predictions differ by about 30 % due to different 
assumptions on the cosmic ray spectrum and composition and the description 
of the hadronic interaction. Absolute atmospheric neutrino spectra in the in¬ 
terval 320 MeV < < 30 GeV for Ue and 250 MeV < E^^ < 10 TeV are 

measured by the Frejus-experiment [167]. The observed neutrino event types 
can be divided by their experimental separation into contained, stopped and 
throughgoing events (Fig.13). For neutrino oscillation searches it is convenient 
to use the ratio /i/e or even the double ratio R of experimental values versus 
Monte Garlo prediction 


n _ (h/^)exp 


(76) 


where /i denotes muon-like and e electron-like events. Here a large number of 
systematic effects cancel out. The above mentioned calculations agree for R 
within 5 % for Ey between 400 MeV and 1 GeV but show a signihcant dif¬ 
ference in normalisation and spectral shape. This effect can mainly be traced 
back to different assumptions on the production of low energy pions from 10 
- 30 GeV p-Air interactions. This might be improved by the results of the 
recent NA56 measurements [142]. Furthermore the predictions can be cross¬ 
checked with atmospheric muon flux measurements which are closely related 
[168-170]. 

The purest sample to investigate are the contained events corresponding to 
Ey ^ 1 GeV. The events are basically due to quasielastic GG and single pion 
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Fig. 14. Momentum spectra of Super-Kamiokande sub-GeV sample a) for the e-like 
events b) for the //-like events and c) the obtained i?-ratio. R stays essentially flat 
over all five bins (from [172]). 





Fig. 15. Observed R-values by different experiments . No oscillation corresponds to 
R=l. 

NC reactions [171]. Unfortunately the relevant cross sections for these pro¬ 
cesses have a relatively large uncertainty in the energy range of interest. By 
far the highest statistics for the sub-GeV region < 1.33 GeV, where 

Eyis is the energy of an electromagnetic shower producing a certain amount 
of Gerenkov-light) is given by Super-Kamiokande . With a significance of 33 
ktxy they accumulated 1158 //-like and 1231 e-like events in their contained 
single ring sample [172]. The capability to distinguish e-like and //-like events 
in water Gerenkov-detectors was verified at KEK [173]. The momentum spec¬ 
tra are shown in Fig. 14. The value obtained with two independent analyses is 
given hj R = 0.61 ± 0.03{stat.) ± 0.05(s//s.). A compilation of experimental 
results is shown in Table 8 (Fig. 15). While Frejus and NUSFX are in agree¬ 
ment with expectations, it can be seen that the water Gerenkov detectors and 
Soudan2 show a signihcant reduction. Besides looking on the R-ratio for os¬ 
cillation searches, the zenith angle distribution can be used (Fig. 16). Because 
the baselines are quite different for downward (L 20 km) and upward going 
muons (L 10^ km), any oscillation effect should show up in a zenith an¬ 
gle dependence. The recent distributions from Super-Kamiokande also for the 
multi-GeV sample {Ey^g > 1.33 GeV), consisting of contained and partially 
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Table 8 

Compilation of existing R ratio measurements. The statistics is now clearly domi- 


nated by Super-Kamiokande . The 

no oscillation case 

corresponds to R = 1. 

Experiment 

R 

stat. significance (kT x y) 

Super-Kamiokande (sub GeV) 

0.63 ±0.03 ±0.05 

33.0 

Super-Kamiokande (multi GeV) 

0.65 ±0.05 ±0.08 

33.0 

Soudan2 

0.61 ±0.15 ±0.05 

3.2 

IMB 

0.54 ±0.05 ±0.11 

7.7 

Kamiokande (sub GeV) 

0.60l[];[]^ ± 0.05 

7.7 

Kamiokande (multi GeV) 

0.57l[];[]^ ± 0.07 

7.7 

Frejus 

1.00 ±0.15 ±0.08 

2.0 

Nusex 

U.yD_Q 28 

0.74 


0.5 

“ -1 -0.5 0 0.5 1 ° -0-5 o 0.5 1 

cos© cose 

Fig. 16. Zenith angle dependence of R for the sub- (left) and multi-GeV (right) 
sample of Super-Kamiokande . The line in the multi-GeV sample corresponds to 
- i^r oscillations with sin?29 =1 and Am^ = 5 • (from [172,174]). 

contained events, are shown in Fig. 17 showing an up-down asymmetry which 
could be explained by neutrino oscillations [174,175]. To verify this assumption 
an L/E analysis for fixed , as the one proposed for the LEP-experiments 
in [176], is done, which shows a characteristic oscillation pattern. From the 
zenith angle distribution and the momentum spectra it seems evident that 
there is a dehcit in muon-like events, which might be explained by — v-r 
or — Vs oscillations. The region allowed by — v^ oscillations is shown in 
Fig.18. An independent three flavour analysis results in a best £t value of Am^ 

~ 7 ■ 10“^eR^ for maximal mixing [177]. Additionally the CHOOZ-result ex¬ 
cludes all Kamiokande data to be due to v^ — z/g oscillations and are shown for 
comparison in Fig. 18 as well. Moreover in a recent analysis of all atmospheric 
data including the earth matter effect (see chapter 6.3.2), the CHOOZ-result 
rules out the z/^ - z/g solution for Super-K at 90 % CL [178]. Furthermore dif¬ 
ferent oscillation channels might be distinguished by a detailed investigation 
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Super-Kamiokande Preliminary 



COS 0 

Fig. 17. Zenith angle distribution of electron- (left) and muon like (right) events for 
the sub- and multi GeV sample of Super-Kamiokande . The bars correspond to MC 
expection without neutrino oscillations and are reflecting the statistical uncertainty. 
The ± 20 % normalisation uncertainty is not shown. The normalisation is adjusted 
to the upward going events. The dotted curve corresponds to a -Vr oscillation 
scenario with si'n?26 = 1 and Am^ = 5 • 10“^eK^. A clear /x-dehcit coming from 
below can be seen. The ratio R obtained from the histograms is shown in Fig. 16 
(from [174]). 

of up-down asymmetries [179] or by measuring the NC/CC ratio [180]. 

Neutrino events at higher energies are detected via their CC reactions pro¬ 
ducing upward going muons . The effective detector area can be increased 
because of the muon range allowing CC in the surrounding rock (see chap¬ 
ter 7.2.2). The corresponding muon flux of the used horizontal and upward 
going muons has to be compared with absolute predictions. One also has to 
take care of the angular dependent acceptance of the detector. Here the main 
uncertainty for the neutrino flux stems from kaon production and the knowl¬ 
edge of the involved structure functions. The behaviour of low-energy cross 
sections is dicussed in [171]. Also here the models can be adjusted to recent 
muon flux measurements in the atmosphere [168] even though one has to take 
into account that for E>100 GeV relatively more neutrinos are produced by 
kaon-decays while the muon-flux is still given dominantly by pion-decay. The 
observations of upward going muons are compiled in Fig. 19. A zenith angle 
distribution from upward going muons as measured with Super-Kamiokande 
is shown in Fig. 19. Two independent ways of verifying the oscillation solu¬ 
tion are the ratio of stopped/throughgoing muons and the shape of the zenith 
angle distribution [181]. Both were done by Super-Kamiokande and support 
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Super-Kamiokande Preliminary 




1 ^ 29 


Fig. 18. a) - Ur allowed region of Super-Kamiokande analysis of zenith angle 
distribution (from [175]). b) u^, - u^ exclusion plot as given by reactor experiments 
especially CHOOZ. This excludes any interpretation of atmospheric neutrino data 
as observed by Kamiokande by oscillations in this channel (from [130]). 


Upward Through Going Muon Flux 




Fig. 19. Zenith angle distribution of upward going muons compiled by Fogli et 
al. [182] as of Aug. 1997 (left). The up-to-date spectrum as observed by Su¬ 
per-Kamiokande is shown on the right. The line corresponds to the flux predictions 
of Honda et al. [183] and using the GRV 94 DIS structure functions [184]. An oscil¬ 
lation scenario with siin?29 = 1 and Am^ = 5 • 10“^eK^ including a factor a = 1.2 
in the normalisation is shown as a dotted curve (from [175]). 


their oscillation evidence [174]. 
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Fig. 20. Energy generation in the sun via the pp und CNO cycles. 

6 Solar neutrinos 


The closest astronomical neutrino source is the sun. The investigation and 
understanding of the sun as a typical main sequence star is of outstanding 
importance for an understanding of stellar evolution. Stars are producing their 
energy via nuclear reactions. The hydrogen burning is done in two ways as 
shown in Fig.20, the pp-chain and the CNO-chain. The net result is the same 
giving 


4p —‘^He + 2z/e + 2Q.73MeV 


(77) 


The prediction of the expected neutrino flux depends on detailed calculations 
of the solar structure resulting in temperature, pressure and density prohles 
and the knowledge of nuclear cross sections for determining the energy gener¬ 
ation. Once the flux is in hand, it is still a matter of detecting this low-energy 
neutrinos typically below 15 MeV with the main component below 500 keV. 
The principle methods are radiochemical detectors using inverse /5-decay and 
real-time experiments looking for neutrino-electron scattering. Because of the 
low cross-sections involved, it is convenient to introduce a new unit for the 
expected event rates in radiochemical detectors called SNU (solar neutrino 
unit) given by 

ISNU = 10“^® captures per target atom per second (78) 


The fundamental equations and ingredients of standard solar models are dis¬ 
cussed hrst. For more detailed reviews see [185-187]. 



6.1 Standard solar models (SSM) 


The sun as a main sequence star is producing its energy by hydrogen fusion 
and its stability is ruled by thermal and hydrodynamic equilibrium. Modelling 
of the sun as well as the prediction of the expected neutrino flux requires the 
basic equations of stellar evolution: 

Mass conservation 


dM{r) 

dr 


47rr^p(r) 


(79) 


where M(r) denotes the mass within a sphere of radius r. 

Hydrostatic equilibrium (gravity is balanced by gas and radiation pressure) 


dp{r) GMir) 

— = — 


(80) 


Energy balance, meaning the observed luminosity L is generated by an energy 
generation rate e 


= 47rrV(r)e (81) 

dr 

Energy transport dominantly by radiation and convection which is given in 
the radiation case by 


dTjr) ^ _3 Kp{r)L{r) 

dr 647r(T 

with a as the Stefan-Boltzman constant and k as absorption coefficient. These 
equations are governed by additional three equations of state for the pressure 
p, the absorption coefficient k and the energy generation rate e: 

P = Pip:T,X), k = k{p,T,X), e = e(p,T,X) (83) 


where X denotes the chemical composition. The Russell-Vogt theorem then 
assures, that for a given M and X an unique equilibrium conhguration will 
evolve, resulting in certain radial pressure, temperature and density prohles. 
Under these assumptions, solar models can be calculated as an evolution¬ 
ary sequence from an initial chemical composition. The boundary conditions 
are that the model has to reproduce the age, luminosity, surface tempera¬ 
ture and mass of the present sun. The two typical adjustable parameters are 
the ‘^He abundance and the relation of the convective mixing length to the 
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pressure scale height. This task has been done by several groups [187-193]. 
Nevertheless there remain sources of uncertainties. Some will be discussed in 
a little more detail. 


6.1.1 Diffusion 

Several experimental evidences strongly suggest a signihcant mixing and gravi¬ 
tational settling of He and the heavier elements in the sun. The long standing 
problem of - depletion in the solar photosphere can be explained if is 
destroyed by nuclear burning processes, which on the other hand requires 
temperatures of about 2.6 TO® K. Such temperatures exist nowhere at the 
base of the convection zone, therefore "^Li has to be brought to the inner 
regions. Also the measured sound speed prohles in the solar interior obtained 
by helioseismological data can be better reproduced by including diffusion 
processes. Therefore these effects were included in the latest solar models. 


6.1.2 Initial composition 

The chemical abundances of the heavier elements (larger than helium) is an 
important ingredient for solar modelling. Their abundances influence the ra¬ 
diative opacity and therefore the temperature prohle within the sun. Under the 
assumption of a homogeneous sun, it is assumed that the element abundance 
in the solar photosphere still corresponds to the initial values. The relative 
abundances of the heavy elements are best determined in certain kind of me¬ 
teorites, the type I carbonaceous chondrites, which can be linked and found 
in good agreement with the photospheric abundances [194,195]. Abundances 
of C,N and O are taken from photospheric values, the abundance cannot 
be measured and is used as an adjustable parameter. 


6.1.3 Opacity and equation of state 


The opacity or Rosseland mean absorption coefficient k is dehned as a har¬ 
monic mean integrated over all frequencies u 


1 

K 


roo 

JO dT 


/c 


oo 

0 dT 


(84) 


where B,y denotes a Planck-spectrum. The implication is that more energy 
is transported at frequencies at which the material is more transparent and 
at which the radiation held is more temperature dependent. The calculation 
of the Rosseland mean requires a knowledge of all involved absorption and 
scattering cross sections of photons on atoms, ions and electrons. The calcu- 
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lation includes bound-bound (absorption), bound-free (photoionization), free- 
free (inverse bremsstrahlung) transitions and Thomson-scattering. Corrections 
for electrostatic interactions of the ions with electrons and for stimulated emis¬ 
sions have to be taken into acount. The number densities rii of the absorbers 
can be extracted from the Boltzmann and Saha equations. The radiative opac¬ 
ity per unit mass can then be expressed as (with the substitution u= hu/kT) 

1 7 15M"^e“/47r"^(e“ — 1)^ , , , 

— =p 7" - r-77 - du 85 

K 7 (1 - e“) Ei cnni + O-sUe 


where denotes the Thomson scattering cross section. 

The most comprehensive compilation of opacities is given by the Livermore 
group (OPAL) [196,197]. It includes data of 19 elements in a wide range of 
temperature, density and composition. The main contribution to the opacity 
in the centre of the sun is given by inverse bremsstrahlung with a few per 
cent contribution of Thomson scattering. A detailed study on opacity effects 
on the solar interior can be found in Tripathy et ah [198]. 

A further ingredient for solar model calculations is the equation of state, mean¬ 
ing the density as a function of p and T or as widely used in the calculations, 
the pressure expressed as a function of density and temperature. Except for 
the solar atmosphere, the gas pressure outweighs the radiation pressure any¬ 
where in the sun. The gas pressure is given by the perfect gas law, where 
the mean molecular weight /i must be determined by the corresponding ele¬ 
ment abundances. The different degrees of ionisation can be determined using 
the Saha equations. An equation of state in the solar interior has to consider 
plasma effects (normally via Debye-Hiickel treatment) and the partial elec¬ 
tron degeneracy deep in the solar core. The latest equation of state is given 
by Rogers et al. [199]. 


6 . 1.4 Nuclear reaction rates 

A detailed prediction of the solar structure and the corresponding neutrino 
flux relies on the nuclear reaction rates [200,201]. Their precise knowledge 
determines the branching in the complex network of reactions and the yields 
of all isotopes. In contrast to typical energies in the solar interior, which are 
in the keV region (Gamow-region), laboratory measurements are normally 
at about MeV and one has to extrapolate down. Because the cross section 
for non-resonant charged particle interactions is steeply falling, it is usually 
parametrized as 


a{E) 


S{E) 

E 


exp{— 27 rri) 


( 86 ) 
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Fig. 21. Cross section for the ^He i^He ,2p)^Ffe at low energies as obtained with 
the LUNA experiment . For the first time data points directly in the Gamow peak 
are experimentally obtained (from [202]). 

where the Sommerfeld-parameter rj is given by ?] = ZiZ 2 e^/Tiv. In cases of no 
resonances, the nuclear or astrophysical S-factor S(E) should show a rather 
smooth behaviour. It is therefore typically the S-factor which is extrapolated 
down to solar energies. Since the energy of the Gamow-Peak is temperature 
dependent, S(E) is for ease of computation expanded in a Taylor series with 
respect to energy 


S{E) = S(0) + S(0)B + + ... 


(87) 


where S'(0),S'(0) etc. are obtained by a fit to the experimental data. A 
compilation of S(0) values for all relevant reactions of the pp-cycle can 
be found in Table 9. Because such extrapolations contain some uncertain¬ 
ties, the idea is to measure the cross section directly in the relevant en¬ 
ergy range. To reach this goal, several additional requirements have to be 
fulhlled, e.g. going underground. In a hrst step this is done by the LUNA 
collaboration building a 50 kV accelerator at Gran Sasso Laboratory to in¬ 
vestigate the ^He i^He , 2 p)^iLe reaction as the hnal step in the ppl-chain 

[ 202 ] . As can be seen in Fig. 21 , the experimental data exceed the theo¬ 

retical expectation of bare nuclei which is due to a still restricted knowl¬ 
edge of screening effects. Further activities will include an upgrade to a 200 
kV and even 2 MV accelerator making the additional measurements of the 
^Be and cross sections possible. New measurements for 

the ^Be (p, 7 )®il cross section at cms energies between 350 and 1400 keV exist 

[203] . Earlier measurements of Kavanagh et al. [204] and Fillipone et ah [205] 
showed a 30 % difference in the absolute value of S(E) in this region. The 
new measurement seems to support the lower S(E) values of [205]. A further 
proposal to investigate this important cross-section exists by using ISOLDE 
at GERN [206]. 
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Table 9 

Compilation of S(0) values relevant for the pp-process. For most reactions the la 
error is given. 


reaction 

pp^ 

d(p,jfHe 
^He {^He ,2p)^Fe 
^He (p,e’'' t'e )'^He 
^He {a,-if Be 
’^Be {v,-ifB 


S(0) (keV -b) 

4.00(1 ± 0.007l|];[|??) • lO-^^ 
0.25 TO-^ 

5400 ± 400 
2.3 -10-20 
0.53 ± 0.05 

U.Uia_Q QQ2 



Fig. 22. Predicted solar neutrino spectrum including thresholds for various running 
and planned experiments (with kind permission of T. Kirsten). 


Table 10 

Flux predictions from four different solar model calculations. 


Flux ) 

BP95 [187] 

RVCD [191] 

CDFLR [193] 

DS96 [192] 

d’uipp) [10^'^] 

5.91 

5.94 

5.99 

6.10 

4)u{pop) [llflcm~‘^s~^ ] 

1.40 

1.38 

1.40 

1.43 

(l)y(f Be) [lO^cm-^s”^ ] 

5.15 

4.80 

4.49 

3.71 

4>yfl'B) [llflcm~‘^s~^ ] 

6.62 

6.33 

5.16 

2.49 

(l>yfl'^N) [lO^cm-^fi-^ ] 

6.18 

5.59 

5.30 

3.82 


5.45 

4.81 

4.50 

3.74 


6.1.5 Neutrino flux predictions 

Several groups are working on the detailed modelling of the sun in order to 
predict accurately the solar neutrino flux and to reproduce the sound speed 
profile measured with helioseismology. A comparison of the different predicted 
neutrino fluxes at the position of the earth is given in Table 10 and Fig.22. 
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6.2 Solar neutrino experiments 


At present results of five neutrino experiments are available, namely the 
chlorine-experiment, the two gallium experiments GALLEX and SAGE 
and the only real-time experiment Kamiokande and its follow-up Super- 
Kamiokande. The discussion follows the historical ordering starting with the 
chlorine-experiment. 


6.2.1 The chlorine-experiment 

The origin of neutrino astrophysics is the chlorine solar neutrino experiment 
by R. Davis in the Homestake mine in South Dakota [207,208]. The detection 
reaction is 


( 88 ) 


with a threshold of 814 keV. Therefore it is basically sensitive to and 
^Be neutrinos with small contributions due to pep, and neutrinos. 
All, except the ^B neutrinos , lead only to the ground state of whereas 
is also populating excited states including the isotopic analogue state. The 
cross section for the reaction (88) averaged over the ^B spectrum has been 
measured recently to be [209,210] 

1.14 ± 0.11 ■ 10"^^cm^ (89) 


The predicted SNU-rate for the experiment due to different flux contributions 
is given in Table 11. The production of one atom/day corresponds to 5.35 
SNU. The experiment consists of 615 t tetrachloroethylene (G 2 GI 4 ) under a 
shielding of about 4000 mwe. The natural abundance of is 24 % resulting 
in 2.2 ■10^'^ target atoms. An extraction of the produced ^^Ar happens roughly 
every two months, and the extraction efficiency is controlled by adding a small 
amount of isotopical pure inert ^^Ar or ^^Ar . To do this, helium is flushed 
through the tank taking the volatile argon out of the solution and allowing 
the collection of the argon in a cooled charcoal trap. After purihcation, the 
argon is hlled with the counting gas PIO into specially developed low-level 
miniaturized proportional counters. The detection reaction uses the EG of 
37 Ar 


e -f ^’^Ar —^'^Cl + Ue 


(90) 


with a half-life of 35 days and focuses on observing the 2.82 keV Auger elec¬ 
trons. To discriminate further against background energy and pulse rise time 
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Table 11 

Contributions (in SNU) of the different flux components of solar neutrinos to the sig¬ 
nal in different radiochemical detector materials of running or planned experiments 
. The fluxes of BP95 are used._ 


Source 

37(7/ 

71 Ga 

127j 

7Li 

i3iXe 

PP 

0 

69.7 

0 

0 

9.7 

pep 

0.22 

3.0 

1.85 

9.17 

1.6 

'^Be 

1.24 

37.7 

13.0 

9.78 

17.8 


7.36 

16.1 

18.4 

25.8 

12.7 

137V 

0.11 

3.8 

0.73 

2.62 

1.6 

150 

0.37 

6.3 

2.43 

13.4 

1.8 

Sum 

9.3 

136.6 

36.5 

60.8 

45.2 



Fig. 23. SNU-rate as function of time for the Cl-experiment (from [211]). 

information are used and the counters are plugged into a special low-level 
shielding. The average measured value using 108 runs after starting the ex¬ 
periment in 1970 is given by [211] 

2.56 ± 0.16(stot) ± 0.15{sys.)SNU (91) 


whereas the single runs can be seen in Fig. 23. The theoretical expectations 
are 9.3 ± 1.4 SNU [187], 4.1 ± 1.2 SNU [192] and 6.4 ± 1.4 SNU [190]. This 
discrepancy is the origin of the solar neutrino problem. 


6.2.2 Real time water Cerenkov detectors 

The only real time solar neutrino experiments are Super-Kamiokande and its 
precessor Kamiokande. The detection principle is the Cerenkov light produced 
in neutrino-electron scattering within the water. Energy and directional infor¬ 
mation are reconstructed from the corresponding number and timing of the 
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Fig. 24. Angular distribution of observed recoil electrons in Super-Kamiokande . 
Clearly visible is an increase in direction of the sun (cos 6 =1). The double angle 
plots (right) with the sun at (0,0) show clearly that the neutrinos are coming from 
the sun (from [213]). 

hit photomultipliers. The scattering angle of the struck electron is related with 
the incident neutrino energy as 

\ _|_ Hie 

(92) 


where Tg denotes the kinetic energy of the recoil electron. The directional 
information is shown in Fig.24. While Kamiokande consisted out of 3000t of 
water using only 680 tons as hducial volume, Super-Kamiokande consists of 
50000 t using 22.5 kt as hducial volume. The detector threshold is 6.5 MeV 
for Super-Kamiokande (in the late stage of Kamiokande it was at 7.5 MeV) 
making these detectors only sensitive to neutrinos. The measured huxes 
are [212] 


<h(®i? ) = 2.80 ± 0.19 ± 0.33 ■ lO^cm ^ Kamiokande (hnal) (93) 
<h(®i? ) = 2.44 ± 0.05(stat.)lQ;Q7(s|/s.) ■ Super-Kamioka4fl4) 

where the theoretical expectations are 6.62 ■ (BP), 4.52 ■ 

10®cm“^s“^ (TCL) and 2.49 ■ 10®cm“^s“^ (DS). The ratio of measured to 
expected electron recoil spectrum is given in Fig. 25. Super-Kamiokande re¬ 
cently implemented a low energy trigger with a threshold of 4.5 MeV. 


6.2.3 The gallium experiments 

The only experiments which are able to measure the pp- neutrinos directly are 
GALLEX and SAGE, using '^^Ga as target material. The underlying reaction 

is 


'^^Ga -I- z/g Ge + e 


(95) 
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Fig. 25. Energy spectrum oi^B neutrinos as observed by Super-Kamiokande divided 
by the solar model of [184]. The solid line corresponds to a fit with a SMA solution 
(from [212]). 

with a threshold of 233 keV. The main difference between the two experiments 
lies in the chemical state of the gallium and therefore also in the extraction 
of the produced '^^Ge . While GALLEX is using 30 t of gallium in form of 
a 110 t GaCls solution, SAGE is using about 57 tons of metallic gallium. 
After extraction, is converted into GeH 4 and filled together with Xe into 
special miniaturised proportional counters. The detection relies on the Auger- 
electrons and X-rays from K and L-capture in the decay producing two 
lines at 10.37 keV and 1.2 keV. As in the chlorine-experiment, besides the 
energy information also pulse rise time analysis is used and the counting is 
done inside a special low-level shielding. Both experiments for the hrst time 
checked their overall efficiency by using MGi ^^Gr sources. The present results 
are [214,215] 

76A±6.3{stat.)tti{sys.) SNU GALLEX (96) 

66.6l®;®(stat.)13;®(s|/s.) SNU SAGE (97) 

with theoretical predictions of 137i'i7 SNU [216], 123 ±8 [190] and 115 ±6 SNU 
[192]. Glearly the experiments are far off. The individual runs for GALLEX 
are shown in Fig. 26. 


6.3 Solutions to the solar neutrino problem 


The observed results split the solar neutrino problem into three. The hrst and 
original one is that the measured rate in the chlorine - experiment, dominated 
by neutrinos , is less than the SSM prediction. This problem might be 
explained as an astrophysical solution by reducing the - hux by a temper¬ 
ature decrease in the solar core or by the involved uncertainties in the nuclear 
cross sections. On the other hand Super-Kamiokande measures the ^B - hux 
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Fig. 26. SNU rate as a function of time for the GALLEX experiment (from [215]). 

and taking this valne as a fact, the contribntion of the nentrinos to the 
chlorine signal already exceeds the experimental valne. This is independent 
of any solar model. There is no astrophysical scenario which could distort 
the beta spectrum in such a way that both experiments are in agreement. 
Any possible deviation is shown to be at maximum at the 10“® level [217]. 
There is almost no room for the ^Be neutrinos , but ^B is produced from 
^Be . The third problem is, that the gallium experiments do not allow any 
signihcant contribution beside the expected pp-value. Also here there is al¬ 
most no room for the "^Be neutrinos contributing more than 30 SNU in the 
SSM. A fit to all available data, independent of a solar model, is given by Hata 
et al. [218] and shown in Fig. 27. The best £t values achieved for the fluxes 
are ^(^Be )/^(^Be )ssm = —0.6 ± 0.4 and )/<h(®i? )ssm = 0.4 ± 0.5. 
Restricting the fluxes to physical regions (<F > 0) changes the result to 
<h(^Re )/^{jBe )ssm < 0.1 and )/<F(®i? )ssm = 0.38 ± 0.05 using the 
solar model of Bahcall and Pinnsoneault [187]. 

6.3.1 Astrophysical and nuclear solutions 

Typical astrophysical solutions try in some way to reduce the central tem¬ 
perature in the sun to account for less high energetic neutrinos . For an 
overview see [186]. Because of the strong temperature dependence of the ^B - 
flux (oc T^®), a reduction to 96 % of the SSM of BP value could explain the 
Super-Kamiokande data. However the ratio <F(^Re )/<l>(®R ) oc increases 
in contrast to experiments which basically imply no '^Be neutrinos at all. This 
is the main reason why neutrino solutions are preferred. A way to circumvent 
this naive T dependence is given by Gumming et al. [220] assuming a slow 
mixing of the solar core on timescales characteristic of ^He equilibrium. The 
result is a remarkably different out-of equilibrium ^He - profile in the solar 
core leading to two consequences: First of all, more helium is produced via 
the ^He i^He ,2p)^Fre chain reducing the Be and ^B neutrino flux and sec¬ 
ondly, the short-living ^ Be is produced at higher temperature favouring the 
Be (p, 7 )®R reaction with respect to Be Li . The combined effect is a 


52 







I 



0.5 


0.0 


• Monte Carlo SSMs 

• TL SSM 
O Low Z 

□ Low opacity 
■ WIMP 
▲ Large S ,, 

O Small S,., 

V Large S 3 , 

♦ Mixing models 

A Dtir-Shaviv model 

Cumming-Haxton model 



Combined data 

w/o luminosity const. . ■' 

95, 99% C.L. , ' ^ 1 


0.0 


0.5 

(t)(^Be)/ (t)(7Be)ssM 


1.0 


Fig. 27. Constraints from the Cl,Ga and Cerenkov-experiments in comparison with 
several flux predictions. In the upper right corner are the solar models from Bahcall 
and Pinnsoneault [187], the dashed curve corresponds to a power law behaviour 
of the solar core temperature and the other mentioned solar models are explained 
in Hata and Langacker [219]. Reducing the Sn factor leads to an evolution in the 
direction of the arrow. The fit to the experimental data results in negative values for 
the ^Be flux (grey area). On the right side no luminosity constraint is implemented 
(from [218]). 



Fig. 28. Isothermal sound speed {u = p/p) profile compared with the reference 
model of [222] and the standard solar model of [187]. Shown is the deviation of 
sound speed Au from the calculated value (from [222]). 

somewhat reduced flux and a significantly reduced ^Be flux. 

All the models experience signihcant constraints from helioseismological data. 
The agreement between the measured and calculated sound speed prohles are 
in good agreement even in the solar interior (Fig. 28) [221]. In the region 0.2 
< t/Rq < 0.65 the deviation from expectation is less than 0.5 % and even in 
the solar core it seems to be less than 4 %. 


6.3.2 Neutrino oscillations in matter 

An elegant solution to the solar neutrino problem are neutrino oscillations 
either in vacuum or in matter (MSW-effect). The latter offers the chance to 
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Fig. 29. MSW-solutions. The excluded region is given by the non-observation of any 
day-night effect. The two islands (grey) are the allowed parameters in agreement 
with all experimental data (from [218]). Also shown are the allowed contours for 
the single experiments. 

suppress the neutrinos of intermediate energies completely, but leaving the 
low-energy neutrinos untouched and the high energy neutrinos only partly 
suppressed. While in solar matter z/g can interact with electrons via charged 
and neutral currents, other neutrino flavours only have neutral current inter¬ 
actions. This leads to different forward scattering amplitudes making the mass 
eigenstates within the sun a function of the electron density. For the oscillation 
amplitude, a resonance behaviour occurs, allowing maximal mixing even if the 
vacuum mixing angle is small [223,224]. The resonance occurs at an electron 
density W of 

/\m?cos29 

Ne = -^- 

2y/2GFE 



In the adiabatic limit, where the electron density along the trajectory of the 
neutrino changes slowly enough, a full conversion is achieved, while in the 
non-adiabatic limit there is a certain transition probability. In case of a lin¬ 
ear change in electron density the probability is given by the Landau-Zener 
probability 


P{E) = exp{ 


'KAm‘^sin‘^ 26 Rs 

E 


) 


(99) 


with Rs = 6.6 ■ 10® m. Because the running experiments have different energy 
thresholds, different contours in the Am® - sin‘^26 plane arise, which overlap 
only in small regions. Careful statistical analyses have been done by several 
authors [218,225-227]. The preferred solutions are a non-adiabatic or small 
angle solution (SMA) at Am® ~ 5 ■ 10“®el7® and sin^26 ^ 0.008 and a large 
angle solution (LMA) at Am® 1.6 ■ 10“®el7® and sm®26* ~ 0.6 (Fig. 29). 
Also vacuum oscillations (VAC) are not ruled out [228], giving parameters 
of Am® 6 ■ 10“^^eV® and sirP29 k, 0.9 (Fig. 30). A hybrid solution of 
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Fig. 30. Am? versus siv?29 plot if vacuum oscillations are the explanation for the 
solar neutrino deficit (from [230]). 

vacuum oscillations and MSW-effect seems possible in a three flavour scenario 
[229]. The upcoming SNO detector will be able to distinguish between all 
the different solutions. By looking for the electron recoil spectrum as well 
as the ratio of charged currrent/ neutral current reactions, they should in 
principle be able to distinguish the solutions with a high conhdence level. Also 
the upcoming radiochemical detectors with their different energy thresholds 
provide further rejection power between the different solutions. The earth can 
regenerate some of the neutrinos for certain parameter values. By looking for 
the resulting day(D)-night(N) effect, Super-Kamiokande gives a value of 

X 100 = -2.3 ± 2.0(stat.) ± 1.4(s|/s.) (100) 

showing no hints for such an effect and excluding certain parameter regions in 
Fig. 29. Because of the long exposure time of the radiochemical experiments, 
they are not able to search for this effect. Including a second flux independent 
quantity A(T)/(T), the deviation of the average measured kinetic energy of 
the electrons from the standard value, Fogli et ah [230] obtain from the binned 
spectrum of Super-Kamiokande 

A(T)/(T) X 100 = 0.95 ± 0.73 (la total) (101) 

resulting in a small preference for the SMA-solution (Fig.31). 

6.3.3 Neutrino magnetic moments in matter 

A similar resonance behaviour can also occur if the solar neutrino problem 
is solved by a neutrino magnetic moment. The pure spin-flavour precession 
^eR ln the solar magnetic held cannot explain the data because it results 
in an energy independent suppression. By allowing spin-havour transitions like 
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Fig. 31. Day-night effect versus the deviation of the electron recoil spectrum. By an 
accurate measurement of these two flux-independent quantities it seems possible to 
distinguish the large (L) and small (S) MSW solution. The no oscillation scenario 
is marked as a star (from [230]). 



Fig. 32. Two resonances can occur for neutrinos coming from the solar interior. 
First they experience the RSFP resonance and afterwards, if the conversion is not 
complete, they go through the MSW resonance. 

i/gL —>■ it has been shown that a resonance behavionr can occnr [231]. The 

transition probability can be written as [232] 


P{h'eL r) = 


(2/i5j 


- y/2GFNeffY + 


sin 


-v/Zlr)(102) 


where D is the denominator of the pre-sine factor and Neff is given by Ne — 
Nn/2 (Dirac) and Ne — N^ (Majorana) respectively. The most general case is 
the occnrence of two resonances (Fig. 32). By transversing the sun, neutrinos 
hrst undergo the spin-flavour precession and afterwards the MSW-resonance. 
Depending on the involved Am^ and E, the predicted conversion probability 
can be quite complicated and detailed predictions for the experiments depend 
on the chosen parameters. It is interesting to note, that in case of adiabacity 
in the resonant spin-flavour precession scenario a MSW resonance never will 
occur. Assuming a maximal magnetic held within the sun between 25 -50 kG 
and a magnetic moment of /ij, = 10“^^ fiB the observed data can be explained 
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if Am^ is within a region of 4 ■ 10 ® — 2 ■ 10 In case both mechanisms 

are at work the Am^ - region is shifted to 10 — 10 and sin‘^26 < 

0.25 by allowing a maximal magnetic held between 15-30 kG. Support for this 
scenario could come from the detection of solar z/g which can be produced via 
^eL —^ ^ i^eR- A detailed discussion can be found in [232], 


6.4 Future solar neutrino experiments 


From the discussion in the last section it seems obvious that the investigation 
of the ^Be region needs special attention. At present several experiments are 
under construction of which some should be able to investigate this region as 
well as hnally solve the problem. 


6.4-1 Sudbury Neutrino Observatory 

The next real-time solar neutrino experiment which will be online is the Sud¬ 
bury Neutrino Observatory (SNO). This detector will use 1000 t of D 2 O and 
is installed in the Creighton mine near Sudbury, Ontario. The big advantage 
of this experiment is a model independent test of the oscillation hypothesis by 
using weak neutral and charged currents. The detection reactions are 


Ue + d^e~+p + p 

(103) 

v + d^i' + p + n 

(104) 

n + e ^ n + e 

(105) 


While the hrst reaction is flavour sensitive, the second is not. To detect the 
neutron in the second reaction, two strategies are envisaged: Cl will be added 
to the heavy water, to use the process and/or a set of He- 

hlled proportional counters will be deployed. The threshold will be around 5-6 
MeV and start of measuring is expected in 1998. Beside the NC/CC ratio, 
the measured electron recoil spectrum of the neutrinos will provide strong 
discrimination power among the different scenarios (Fig.33). 


6 . 4.2 BOREXINO 

An experiment especially designed to explore the intermediate region contain¬ 
ing the Be line at 862 keV is BOREXINO. It will use 300 t of scintillator. 
The detection relies on neutrino-electron scattering and the detector will be 
sensitive to neutrino energies larger than 450 keV. Of special importance is 
the produced ” Compton-Plateau” from the 862 keV line with an expected 
event rate of 50 per day according to the SSM or about 10 per day in case of 
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Fig. 33. Statistical power of Super-Kamiokande and SNO to distinguish between 
the different solutions of the solar neutrino problem. Shown are the moments of the 
electron sprectrum (for details see [225]). The SMA can be checked at the 3a level, 
the error bars correspond to the experimentally allowed MSW-region (with kind 
permission of J.N. Bahcall). 

the MSW-effect. The detector is extremely sensitive to impurities requiring a 
background level of 10“^® g(U,Th)/g. The ability to achieve such low back¬ 
ground levels could be demonstrated in a smaller pilot experiment (CTF). 
BOREXINO is currently installed in the Gran Sasso Laboratory. Data taking 
should start around 2001. 


6 . 4.3 ICARUS 

The ICARUS experiment plans to use in full scale about 3000t LAr in form 
of a TPC for solar neutrino detection. The technique offers precise tracking 
with high resolution dE/dx measurement, full event reconstruction in a bubble 
chamber-like quality and full sampling electromagnetic and hadronic calorime¬ 
try. Besides neutrino-electron scattering with a possible threshold of about 5 
MeV, also the capture to the isotopic analogue state in 

z/e + ^°Ar ^e-+^^K* (106) 


will be used. The threshold for this reaction will be 5.9 MeV and allows the 
detection of photons with a sum energy of 4.38 MeV from the de-excitation of 
in coincidence with the electron. By comparing both reaction rates, it is 
possible to get direct informations on the oscillation hypothesis. A hrst 600t 
module will be installed in Gran Sasso Laboratory 1999. 


6 . 4.4 HELL A Z 

A large real time experiment even for the detection of pp-neutrinos is the 
proposed HELLAZ experiment. It will consist of a 12 t high pressure helium 
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TPC at liquid nitrogen temperature. A smaller prototype for demonstration 
of the basic principle exists. 


6.4.5 127 / 

A new detector relying on jg at present installed in the Homestake mine 
near the chlorine experiment [211]. The detection reaction is 

z/e +127 I ^e- +127 Xe (107) 


with a threshold of 789 keV. All the experimental mechanisms like extrac¬ 
tion and detection are similar to the chlorine experiment . The extraction will 
happen in two cycles (day/night) to two different charcoal traps with an ac¬ 
cumulation time of about 1 month. i27Xe will decay to excited states of i27/ 
making a coincidence detection of the gamma with the Auger electrons possi¬ 
ble. The detector will use 235 t of Nal solution with a total of 100 t iodine and 
will be in operation soon. The expected event rate is about 36.4 SNU where 
14 SNU results from Be neutrinos . In contrast to the chlorine experiment , 
there is no bound isotopic analogue state in the I-Xe system, increasing the 
relative contribution from ^Be with respect to ^B . 


6 . 4.6 Li 

In comparison with chlorine it might be interesting to envisage a 
Li experiment. Both have nearly identical thresholds but the contribution 
of 7_Be and ^B neutrinos to the signal is qnite different. Snch a detector is 
using the reaction 


Li + z/g —> Be + e 


(108) 


and is currently under construction [233]. The energy threshold is 860 keV 
making this detector sensitive to ^Be and ^B neutrinos. The plan is to use 10 
tons of metallic Li. The solubility of Be in Li decreases with falling temper- 
atnre, making an extraction with cooled traps possible. Only 10.4 % of EC 
from ^Be prodnce a 478 keV photon which could be used for detection by 
conventional techniques. Therefore cryogenic detectors are necessary to mea¬ 
sure the Auger electrons and nnclear recoil adding to an energy of 112 eV. A 
prototype of 300 kg Li is constructed and presently under investigation. The 
predicted rate is 60.8lg SNU [233]. 
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A different approach to use low threshold detectors allowing solar neutrino 
spectroscopy is given by Raghavan [234], The idea is to use an isotope I and 
populate excited states in the neighbour isotope F which will decay after a 
short time (10 - 100 ns) with y-emission (z/e+/ —> e~+F* F+'j). This makes 
the use of delayed coincidence techniques possible. To prevent the mother 
isotope from decaying by single beta decay, the idea is to use double beta 
decay candidates like or . By using different excited states, 

it is even possible to compare different contributions of the solar neutrino flux. 
A 15 % Yb loaded scintillator could be successfully created, still fulhlling all 
experimental requirements. A lOOt scintillation detector containing lOt Yb is 
taken into consideration. 


6 . 4.8 

Another idea is to use the reaction 

+ z/g ^ ^ g- (1Q9) 


The threshold would be 352 keV and the expected rate is about 45 SNU 
[235]. For the different contributions see tab. 6.2.1. A 1 kton detector would 
result in a detection rate of 1500 events/year according to the SSM, where 
^Be neutrinos would contribute 37 % of the signal. A liquid xenon detector 
like those proposed for dark matter searches would be the appropriate choice. 


6 . 4.9 GNO 

Because GALLEX is hnished and the wish to continue with measurements 
of pp-neutrinos over the next decade, it was decided to continue as Gallium 
Neutrino Observatory (GNO). For the hrst 2 years it will continue with the 
same setup as GALLEX but then an increase of the gallium mass to 60 t and 
even 100 t as well as technical improvements are foreseen. 


7 Astrophysical aspects of neutrinos 


7 .1 Neutrinos from supernovae 


Among the most violent stellar events are supernova explosions. Supernovae 
are products of the late phase of stellar evolution and can be divided in two 
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classes. Supernova Type la are C,0 white dwarfs in binary systems, which 
exceed the Chandrasekhar-mass 

Mch = 5.72 Y^Mq ( 110 ) 

by accreting matter from a companion, leading to a thermal deflagration of the 
white dwarf. With this type of supernova no neutrino emission and hydrogen¬ 
lines are connected. Stars with masses of M ~ 8Mq burn nuclear fuel up to 
iron group elements. If the iron core passes the Chandrasekhar-mass, it will 
collapse to a neutron star or even black hole, creating a supernova explosion. 
Because of the ejected outer hydrogen shell, hydrogen-lines can be observed 
and the supernova is called Type II. For a more detailed classihcation scheme 
see [236]. Detailed discussions on the mechanism of supernova explosions can 
be found in [237,238], only the principal scheme relevant for neutrino physics 
is outlined here. Typical values at the beginning of the collapse are a central 
density of p ~ 4 ■ 10^gcm~^, a temperature of 8 ■ 10® K and an electron per 
nucleon fraction of W ~ 0.42. The gravitational force is basically balanced by 
the pressure of degenerated electrons. Photo-disintegration of iron via 

13^He + An-12AA MeV (111) 

and electron captnre on free protons and on heavy nnclei 

e~+p^n + h'e, e~A A + Ue ( 112 ) 

rednce the electron density. Therefore the star loses its pressnre snpport and 
it collapses. This collapse stops when the iron core reaches nnclear density 
because of the now stiff eqnation of state . Because of an overshoot this part 
of the core bonnces back with an energy depending on the nnknown eqnation 
of state for overdense nnclear matter. The onter part of the core still continnes 
to fall in, thus producing a pressure discontinuity at the sonic point, which 
develops into an ontgoing shock wave. Depending on the energy of the shock 
wave, it is able to escape the iron core and to create the explosion (’’direct 
explosion mechanism”) or it stalls in the core and needs some additional energy 
input to be successful (’’delayed explosion mechanism”) (Fig.34). The total 
binding energy released in such an event is of the order 5 ■ 10^^ erg, where 
nentrinos carry away about 99 %. 


7.1.1 General remarks 

The emitted neutrino spectrnm consists basically of two parts. The first one 
is a peak of Ue resnlting from the deleptonisation period (eq. 112) and lasting 
only a few ms. Dnring the collapse phase the core becomes opaqne even for 
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Fig. 34. Time development of a type II supernova explosion after core bounce. As the 
nuclear matter is over-compressed in the collapse, a rebounce occurs and produces 
a shock wave (from [239]). 

neutrinos and they diffuse within the core. The outgoing shock wave dissoci¬ 
ates the infalling iron nuclei, increasing the mean free path for the neutrinos 
which pile up behind the shock. When the shock traverses the neutrinosphere 
(dehned in a similar sense like the photosphere of the sun) all this v^. will be 
emitted at once. The mean energy is relatively high (Tli,) ^ ISMeV, but the 
total energy released is only about 10^^ erg. The second contribution comes 
from the Kelvin-Helmholtz cooling phase of the proto-neutron star acting as 
blackbody source and producing neutrinos dominantly by bremsstrahlung. All 
flavours are emitted in more or less equal numbers (Fig.35) within a few sec¬ 
onds. Because and have no charged current interactions they have a lower 
opacity and decouple at higher temperature and density. Also the opacity for 
z/e is lower than for z/g because less protons are available and the opacity is 
dominated by z/g + n —p-|-e“ and Ue + P ^ n + e'^ respectively. Therefore one 
typically hnds {Ey) = 10-12 MeV (z/e ), {Ey) = 14-17 MeV (z/e ) and {E^) 
= 24- 27 MeV (z/^ ,z/.r )• Because the energy is approximately equipartitioned 
between the flavours, the fluxes behave as ^(z/e) > *h(z/e) > <F(z/^, z/.^)- The 
neutrino spectrum might well be described by a Fermi-Dirac shape 

dLp ^ Ey ( 113 ) 

dEy 1 -k eXp{Ey/Ty - Tjy) 

where rjy is a degeneracy parameter. A detailed study of the time evolution 
of the neutrino luminosity is rather complex and depends on many parame¬ 
ters like the equation of state , the mass of the collapsed core at bounce, the 
amount of postbounce accretion and the temperature prohle after collapse. So 
far the only event which could be detected is SN 1987A . 

In addition past supernovae could create a relic neutrino background which 
is discussed in [241,242]. The total integrated flux predicted by Malaney[241] 
using the redshift evolution of Hl-gas and the associated star formation rate is 
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Fig. 35. Time evolution of the neutrino luminosity and average energy of a supernova 
explosion model from [240]. The core bounce is 3-4 ms before the neutronisation 
burst of Vf. (from [240]). 

between 2-5.4 which is a factor of 10 less than previous estimates of 

[243,244]. The predicted spectrum has a peak in the region 2-5 MeV which is 
signihcantly due to supernovae with redshifts larger than 1. Unfortunately the 
background of terrestrial reactors, solar and atmospheric neutrinos will make 
a detection in this region very unlikely. Experiments like Super-Kamiokande 
therefore have to search in the region 15-40 MeV, where the theoretical pre¬ 
dictions for the supernova background fluxes are approximately the same. 


7.1.2 SN 1987A 

On February 23, 1987 the blue supergiant Sanduleak -69 202 exploded in the 
Large Magellanic Cloud at a distance of about 50 kpc [245]. For the hrst time 
neutrinos could be detected, bringing theoretical supernova model calculations 
in the regime of experimental verihcation. Four detectors claim observation of 
the neutrino burst, namely IMB [246,247], Kamiokande [248,249], the Baksan 
Scintillator Telescope [250] and the Mont Blanc Liquid Scintillator detector 
[251]. The Mont Blanc detection happened about 5 h earlier than the detec¬ 
tion of the other experiments. Because of the relative low energy of the events, 
the non-observation of any signal during this period in the much larger water- 
Cerenkov detectors and a missing astrophysical scenario for producing two 
neutrino bursts, this detection is normally considered as a background fluctu¬ 
ation. 

The relevant interaction processes in the water detectors are VeP ne~^, 
ve — ve elastic scattering and -t- e~ which becomes the dom¬ 

inant contribution for z/g for Ey larger than about 30 MeV (Fig.36). For the 
scintillator detectors the reaction is absent but at Ey larger than about 
30 MeV the reaction ^ -|- e~ contributes. By far the largest cross 

section is UeP —> ne~^ resulting in an isotropic event distribution and is sug¬ 
gesting that all observed events are due to interactions. 

The observed numbers of neutrinos are 11 events within 12 s (Kamiokande), 8 
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Fig. 36. Different cross sections involved in the detection of supernova neutrinos. 
By far the largest cross section is i>eP ne'^ (from [252]). 

events in 5.5 s (IMB) and 5 events in 14 s (Baksan). Some events were already 
attributed to background and are not included. A recent detailed maximum 
likelihood analysis was done by Loredo and Lamb as described in [252]. Their 
best-ht values are 16.9 events plus 5.6 background for Kamiokande, 4 events at 
IMB and 1.8 plus 1 background at Baksan. A two component cooling scheme 
consisting of the Kelvin Helmholtz cooling plus a low energy component which 
mimics the neutrinos emitted during the stalled shock-phase results in a neu- 
trinosphere of 18 km and a total binding energy of 3.08 ■ 10^^ erg in good 
agreement with theoretical expectations. 

The observed signals contain some ’’anomalies”. First of all there is a discrep¬ 
ancy in the neutrino energies between Kamiokande and IMB which imply a 
harder spectrum for IMB. Because of the rather high threshold, IMB is sen¬ 
sitive to the high-energy tail of the assumed neutrino spectrum, which might 
have substantial uncertainties. A second point is the large 7s gap between the 
hrst 8 and last 3 events of Kamiokande. This might be a statistical fluctuation 
because IMB and Baksan do have events in this period. The most striking 
is the deviation from isotropy which is expected if all events are due to Ue 
interactions (Fig.37). The observed distribution is only at the per cent level 
in agreement with isotropy and an explanation for this fact is still missing. 
Concerning neutrino properties several things could be learned even by the 
low statistics of the observed event numbers. Because of the observed pulse 
length, mass effects on the propagation from SN 1987A to the earth restrict the 
neutrino mass to a conservative upper limit of about 25 eV, which could be im¬ 
proved under certain model dependent assumptions down to 13 eV. Because 
the measured number of events corresponds to the theoretical expectation, 
rather stringent bounds on neutrino decay can be given, implying 

—r,, > 5 ■ (114) 

rriy 

The radiative decay channel is additionally limited by the Solar Maximum 
Mission which did not observe any signal which could be related to the neu- 
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Fig. 37. Energy and observed angle as a function of time for IMB (white), 
Kamiokande (black) and Baksan (stars). 

trino burst of SN 1987A [253]. This was confirmed by a systematic search 
for the radiative decay of neutrinos coming from Type II supernova using 
COMPTEL data [254]. A stringent bound also exists on the electric charge of 
the neutrinos. An electric charge would cause longer travel distances for lower 
energy neutrinos because effects due to the galactic magnetic held are more 
important to them. Therefore the charge of the neutrinos can be bounded by 
Q, < 10-i"e. 

Neutrino oscillations and the MSW-effect might also play a signihcant 
role in supernovae. Oscillations in the channel z/g - vx are able to re¬ 
duce the prompt signal signihcantly. A detailed analysis, assuming Am^ 
~ 3 ■ {Ey/IQMeV) to allow for a resonance outside the neutrinosphere 

, reveals that a probability for a conversion of more than 50 % already occurs 
for Am‘^sin^26 ~ 4 ■ eV"^ /IQMeV) extending to large Arn? values. 

Furthermore if the resonance lies outside the neutrinosphere and within the 
shock wave (which is valid for neutrino masses in the region of order 10 eV 
or above and therefore a cosmological interesting region) the higher energy 
and Vr 1 if converted to z/g could help to revive the stalled shock [255]. The 
minimum mixing angle necessary is sin‘^29 ~ 10“® E^y /lO MeV. This would 
imply a reduction of the prompt z/g pulse. The oscillation of z/g could further¬ 
more influence the creation of r-process events. The hot bubble between the 
settled neutron star and the escaping shock wave a few seconds after the core 
bounce seems to be a reliable place for r-process nucleosynthesis. The p/n- 
ratio in the bubble is governed by the reactions z/gU pe~ and ne~^. 

Because neutrinos are much more common than electrons and positrons, the 
p/n ratio is ruled by the spectra and fluxes of the neutrinos . This normally 
results in a neutron rich environment, because z/g are more energetic than z/g 
. Oscillations of z/^ and z/,- to z/g outside the neutrinosphere could make the z/g 
flux more energetic than the z/g . Even a 10 % oscillation effect would drive the 
medium to a proton-rich state. The parameters obtained to allow the r-process 
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also lie in the cosmological region requiring ~ 2 eV and sin‘^26 ~ 10“^ 
for Ve - i^T oscillations [256]. At smaller Am^ the oscillation effect has no 
impact on the r-process because it occurs at too large radii. 


7.1.3 Experimental status 

The experimental observation of SN 1987A launched several new searches 
for supernova neutrinos. Besides specially developed detectors, basically all 
new real-time solar neutrino detectors like Super-Kamiokande , ICARUS and 
SNO will be able to see such neutrinos. The predicted count rate for Super- 
Kamiokande of a galactic supernova at a distance of 10 kpc is 5000 - 10000 Pe 
interactions! Such an event would open the possibility to explore the mass of 
Ufj, and Ur down to 50 eV by using the neutral current excitation of [257]. 
SNO might be especially valuable because of the NC desintegration of D (eq. 
104) which will be dominated by P^ and Ur ■ Detectors like LVD and MACRO 
also offer additional discovery potential. The main component of all detectors 
will still be the Pg detection. A completely different scheme which is mainly 
sensitive to and Ur is the SNBO idea [258] recently put into a realistic de¬ 
tector design in form of OMNIS [259]. The basic idea of [258] is NC excitation 
of nuclei 


uA{Z, N) A*{Z, N) A{Z, N - 1) + u + n (115) 


via the de-excitation by neutron emission. As target material large under¬ 
ground areas of rock or salt should be equipped with neutron detectors. The 
efficiency can be increased by using caverns in the rock for neutron detection. 
With about 200 tons of a Gd loaded scintillator, event rates of more than 2000 
neutrino interactions for a galactic supernova in a distance of 8 kpc seem fea¬ 
sible. Because of their higher energy, the signal would be dominated by and 
Ur interactions. An extension to extragalactic supernovae (an increase to about 
4 Mpc in sensitivity would imply about 1 supernova per year) unfortunately 
seems unrealistic at present times. 


7.2 Neutrinos from other astrophysical sourees 


After describing experimentally observed astrophysical neutrino sources like 
the sun and supernovae there might be other sources of neutrinos of even 
higher energies {E^, > lOOMeU). Their existence may be closely related to 
the recently discovered TeV- y-sources and the sources of cosmic rays. Neutri¬ 
nos are produced via pp-collisions or photoproduction in cosmic beam-dump 
experiments due to the decay of the created charged pions and kaons. The as¬ 
sociated production of neutral pions allows a relation between expected photon 
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and neutrino fluxes (eq.118). The threshold for pion-photoproduction is rather 
high and the minimal proton energy required is given by 


E, = = 7 ■ 10^'>(E,(eV))-^(eV) (116) 


Using the cosmic microwave background (CMB) as target photons, a threshold 
energy for the proton of ~ 5 ■ 10^® eV results (GZK-cutoff) [260,261]. In 
contrast to photons and protons, neutrino propagation is not influenced by 
attenuation or deflection by magnetic fields. Neutrinos give direct information 
on the source location and might reveal sources which have no y-counterpart. 
Even more, while high energy photons are influenced by 77 —> e’''e“ reactions 
on the CMB and in the TeV-region by the same reaction with the presently 
unknown IR-background, they have a limited range for detection, whereas 
neutrinos can be observed to largest distances. 

The main detection reaction on earth will be upward going muons discussed 
in more detail in chapter 7.2.2. 


7.2.1 Sources and predictions 

Several galactic and extragalactic sources are discussed for highest energy 
neutrinos most of them are also investigated for creating and accelerating 
cosmic rays . For a detailed discussion of sources see [262]. The spectral shape 
of the primary cosmic rays follows a power law according to 

$(E) oc (117) 


where the spectrum observed on earth is characterised by 7 ~ 1.7 (for E 
< 10^® eV) then steepening to 7 ~ 2 (’’knee region”) and then changing 
again at 10^® eV (’’ankle”). This spectrum is steeper than the accelerated 
one because of the energy dependence of cosmic ray diffusion in the galaxy. 
Neutrinos produced in interactions with the interstellar gas should follow the 
shape of the primary cosmic ray spectrum up to highest energies. On the other 
hand, if the production occurs at the acceleration site, there is no diffusion 
effect and the spectral shape follows the hard source spectrum (7^2 — 2.2). 
From the above mentioned interaction mechanisms it is clear that there are 
point sources in the sky and a diffuse component. Two examples of possible 
galactic point sources are X-ray binaries (a compact object like a neutron 
star or black hole is accreting matter from a non-compact companion) like 
Cygnus X-3 or young supernova remnants. To produce a detectable signal of 
a few upward going muons per year in a lO^m^ detector. X-ray binaries have 
to convert almost all energy in the acceleration of protons. The supernova 
remnants on the other hand are rare events and typically produce a signal 
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only during a period on the scale of years after the explosion. A guaranteed 
source for neutrino production is the galactic disc, where a diffuse photon 
background due to interactions of cosmic ray protons with interstellar matter 
could be observed by EGRET on CGRO [263]. If neutrinos are coming from 
TT-decay like the observed photons the neutrino flux can be determined by 
[262] 





1 

1-A^ 


( 118 ) 


where is the energy dependent photon attenuation dominated by 77 — 
e’''e“. A second source is our sun, because of cosmic ray interactions within the 
solar atmosphere. Moreover the sun could trap neutralinos y as dark matter 
in its interior. Their yy-annihilation can be a source of high energy neutrinos. 
Predictions for the neutrino -flux on earth for a 500 GeV neutralino y are of 
the order ^ 2 ■ . 

The most prominent extragalactic source candidates are active galactic nuclei 
(AGN) and gamma ray bursters (GRBs). The present -here simplified- pic¬ 
ture of AGN consists of a central supermassive black hole (~ 10® — IO^^Mq) 
accreting matter near its Eddington limit. The accelerated matter will form 
a hot, dense plasma and will be partly sucked into the black hole and partly 
redirected by magnetic fields forming two opposite directed jets perpendicular 
to the accretion disc. In this scenario there are basically two ways to accelerate 
particles by shock acceleration and to produce a neutrino flux. The first possi¬ 
bility is near the central engine as described in [262]. Energy losses take place 
due to processes like pj A+ — mr^ and pj ^ p + + e~ in the radiation 

held as well as pp-collisions in the gas. Both processes give rise to photons 
and neutrinos , but the produced photons from 7r°-decay cascade down and 
are released as X-rays, because the central region is optically thick for energies 
larger than ~ 5 MeV. Grucial for this mechanism to work are the assumptions 
on proton propagation and conhnement in the core region. Note however, that 
the majority of X-ray emitting AGN does not show a nonthermal X-ray spec¬ 
trum (cascade origin) but a thermal spectrum peaking at ~ 100 keV. The 
second source which might explain the observation of TeV-photons from seve¬ 
ral extragalactic sources like Markarian 421 are the highly relativistic jets. It 
is of outstanding interest to know whether the observed photons are created 
from inverse Gompton scattering or synchrotron emission of accelerated elec¬ 
trons or from photoproduction of pious from accelerated protons [264]. The 
observation of neutrinos would help to clarify the situation. 

By integrating over all cosmological AGN, one should also see a diffuse back¬ 
ground of neutrinos in the same way as the diffuse 7 -ray background is ob¬ 
tained. According to Stecker [265] it remains flat up to about 10^ GeV and 
then starts to fall steeply. For energies larger than about 3-10^ GeV it becomes 
dominant with respect to atmospheric neutrinos. While different models agree 




Fig. 38. Predicted diffuse isotropic neutrino flux from the sum over all active galac¬ 
tic nuclei predicted by various models. Also shown is the horizontal and vertical 
atmospheric neutrino flux, dominating in the region below 1 TeV. Two different 
models for prompt neutrinos from charm decay of TIG [160] and model D of Zas et 
al. [271] are also shown (from [270]). 

more or less in their prediction at the high energy end of ~ 10^'^ GeV, orders 
of magnitude differences exist in the lower energy region around 10® GeV. Pre¬ 
dicted fluxes of several models can be found in [266-270]. A combined high 
energy neutrino spectrum for point and diffuse sources can be seen in Fig. 
38. The typical estimates for their detection via upward going muons are in 
the order of 0.1-10 events per year for a 0.1 km^ detector with rather large 
uncertainties. A flux limit of d^/dEy < 7■10~^^GeV~^cm~^s~^ sr~^{90%CL) 
for in the energy region of E^, ^ 2.6 TeV exist from the Frejus-experiment 
[272], already ruling out the model of Szabo et al. [267] and the maximal flux 
predictions of Bhattacharjee et al. [273]. 

Other partially more exotic sources of high energy neutrinos might exist. Neu¬ 
trinos associated with GRBs are discussed in Waxmann et al. [274]. Further 
scenarios are annihilation or decay of superheavy particles like the neutralino 
X [275]. Also evaporating black holes and radiation from topological defects 
like cosmic strings might be sources for neutrinos [273,276]. 


7.2.2 Experimental search 

From the flux estimates of the last chapter, it is rather clear that very large 
detectors are required. The proposed detection reaction for the VHE-neutrinos 
are GG reactions producing upward going muons . The corresponding cross- 
sections are shown in Fig. 39. The GG cross section is given by 

y-T =- n2 I 2 Q ) + Q )(1 - d) ) ( 119 ) 

dxdy TT 
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Fig. 39. Total as well as NC and CC cross sections for high enery neutrinos (left) 
and antineutrinos (right) (from [277]). 

with q{x,Q‘^),q{x,Q‘^) as the quark distribution functions, the Bjorken vari¬ 
able X = Q'^/2mu, m the nucleon mass and y = {Ey — E^)/E^. Towards 
higher energies the contribution from the presently unknown small x-region 
is becoming more important and theoretical extrapolations have to be used 
[277]. To face the background from the atmospheric neutrino flux, which is at 
lower energies typically a factor of 10® larger than cosmic sources (but falls off 
steeper towards higher energies), only muons transversing the detector from 
below can be used (upward going muons). Therefore cosmic sonrces of interest 
have to stay below the horizon for a signihcant time. 

The effective detector size for detecting npward going mnons can be en¬ 
larged because the surrounding material can be used as an additional target. 
The energy loss rate of mnons dne to ionisation and catastrophic losses like 
bremsstrahlnng, pair-prodnction and hadroprodnction is given by 


dEf, 

dX 


E„ 


= —a 


( 120 ) 


with the critical energy e = ^ 500 GeV. This results in a range R of (.^ ~ 

2 km.w.e.) 


R ^ ^ln{l + 



( 121 ) 


For high energy mnons the radiation losses are dominant, resulting in a change 
of energy dependence of the range from linear to logarithmic. For a mnon with 
initial energy larger than 0.5 TeV the range exceeds 1 km. For nentrino ener¬ 
gies larger than abont 40 TeV the interaction length becomes smaller than the 
diameter of the earth resnlting in a shadowing effect. 

The existing nnderground detectors like MACRO, Snper-Kamiokande and 
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LVD have too small sizes to measure statistically significant signals. Therefore 
natural resources of water and ice are used to detect the Cerenkov-light of the 
passing muon. After the termination of DUMAND, the Lake Baikal experi¬ 
ment (water) and AMANDA (ice) are the most advanced projects. 

The Baikal neutrino telescope (NT) is located in a depth of 1.1 km at a dis¬ 
tance of 3.6 km from the shore [278]. The final setup of NT-200 is finished 
recently and consists of 8 strings forming a heptagonal array of seven strings 
around a central string. It consists of 192 pairs of 8 ” photomultipliers arranged 
in a distance of 5 and 7.5 m along the strings. This rather small spacing allows 
a relatively low energy threshold of about 10 GeV. Besides the low threshold 
the main advantages of the experiment are the relatively cheap deployment 
of tubes, because the frozen lake offers a good platform for deployment and 
Lake Baikal is a sweet-water sea containing no whose decay would pro¬ 
duce background. The effective area is between 1000 — 5000m^ depending on 
energy and an increase of the effective area by a factor 20-50 is under consid¬ 
eration. Clear upward going muons have been observed. 

AMANDA [279,280] is at present operating 300 8 ” photomnltipliers in the 
antarctic ice in a depth of 1500 - 2000 m. The basically bubble-free ice offers 
extraordinary optical properties, e.g. the absorption length for A ~ 500 nm is 
around 100 m and the scattering length abont 25 m. The detector offers an 
effective area of lO^m^ with a mean angular resolution of 2.5°. An upgrade 
to a 1 km^ detector (ICECUBE) which would consist of about 5000 photo¬ 
multipliers mounted on 80 strings with a spacing of about 100 m is under 
consideration. Good candidates for upward going muons have been observed 
by AMANDA. 

Two further water experiments in the Mediterranean are in a kind of prepa¬ 
ration phase, namely NESTOR [281] near Greece, in a depth of 3800 m and 
ANTARES [282] near Toulon (France), in a depth of 2000m. 

Associated with two of the projects mentioned above are alternative exper¬ 
iments using different detection techniques. They are called RICE [283] (to¬ 
gether with AMANDA) and SADCO [284] (together with NESTOR). They 
rely on detection of radio and acoustic signals produced ii ^ 1 PeV. In the 
acoustic case, the produced shower particles in the z/g - interaction lose energy 
through ionisation leading to a local heating and density change. The density 
change propagates as sonnd waves through the medium and can be detected 
with an array of detectors like hydrophones allowing also a reconstruction of 
the event by measuring the arrival times and amplitudes. The second method 
nses the effect of coherent radio Cerenkov radiation which is produced as long 
as the wavelength is large with respect to the spatial extension of the shower 
[285]. The nentral pions create an electromagnetic shower of size Im in ice via 
their 7 r° —77 decay, therefore prodncing freqnencies in the region 100 Mhz - 
1 GHz. 

For detection of very high energetic nentrinos the reaction 

z7e + e“ — W~ —> hadrons (122) 
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Fig. 40. Glashow resonance in the z^e cross section. The curves correspond, 
in the low energy region from highest to lowest, to {i)i'e& hadrons 

{ii)Vfj,e iiVe{ni)vee Vee{iv)9ee ^ 9ee{vi)Ufj,e h'f_^e{vii)v^e v^e 

(from [277]). 

can be used which shows a resonance behaviour (Glashow - resonance) at 
s = m^, meaning = 6.3 ■ 10® GeV. The z/gC cross section at the resonance 
is about a factor 30 larger than the corresponding z/N cross section (Fig. 40). 
The held of high energy neutrino astrophysics is still in its initial phase but 
will provide important information in the future. 



7.3 Relic neutrinos 


The thermal history of the universe according to the big bang model pre¬ 
dicts not only a photon background but also a neutrino background. While 
the photons are observed as the cosmic microwave background, the neutrino 
background is still undetected. The temperatures of both are related by the 
relation 


n = (T)iT, (123) 

where T.y is measured quite accurately by GOBE to be = 2.728 ± 0.004 
[286], thus predicting a neutrino background temperature of 1.95 K. The total 
number density and matter density (the havour densities are one third) are 
then given by 


Pu 


8^7 


= -^n.y = 337cm ® 

22 ^ 

(124) 

1 

1 

-3p^ = 0.178^ 

11 cm-^ 

(125) 
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Fig. 41. Neutrino contribution to H as a function of the neutrino mass. Two regions 
remain for stable neutrinos to be cosmological of interest. Either below 100 eV 
(acting as hot dark matter) or above 5 GeV (acting as cold dark matter). Otherwise 
they have to be unstable because they would overdose the universe. 

with the statistical weights = 2, = 2 for Majorana neutrinos and light 

{rriy ~ 300 keV) Dirac neutrinos , otherwise g^ = 4. The mean energy of the 
neutrinos today is 5.28 ■ 10“^ eV making a detection extremely difficult. The 
neutrino contribution to the matter density is given as 

= — = 5.32-10“=^^^ (126) 

Pc eV 


where pc is the critical density and h = Ho/lOOkms ^Mpc ^ the normalised 
Hubble-constant. For rrii, ~ 1 MeV neutrinos become non-relativistic and their 
density is suppressed by a Boltzmann-factor. The behaviour of Qi, as a func¬ 
tion of rriiy is shown in Fig.41. As can be seen, stable neutrinos only exist for 
rriiy ~ 100 eV and for TUiy ~ 2(5) GeV in the Dirac (Majorana ) case (Lee- 
Weinberg-bound). Both allowed mass regions offer massive neutrinos as good 
dark matter candidates, either as hot dark matter (neutrinos in the eV-range) 
or as cold dark matter (neutrinos in the GeV range) [287]. Pure cold dark mat¬ 
ter models predict to many galaxy clusters, which can be avoided by including 
a hot component. Because of their free streaming in the early universe, light 
neutrinos wash out perturbations on small scales, reducing the power there. 
A mixed hot and cold dark matter model with D = 1 and 5 eV neutrinos con¬ 
tributing Qi, ^ 0.2 for h = 0.5 seems to be a good description of the GOBE 
normalised power spectrum and the observed large scale structure. This can 
be improved if the hot component consists of two neutrino flavours having a 
mass of about 2.5 eV [288]. GeV neutrinos as cold dark matter are bounded by 
double beta decay experiments excluding Dirac - neutrinos with standard in¬ 
teractions between 26 GeV and 4.7 TeV as the dominant component [289,290]. 
The linear contribution for light neutrinos to Qi, (eq.l26) can be converted in 
a neutrino mass bound. The condition not to overdose the universe requires 
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for stable neutrinos 


^ 94eVQ^h^ 


(127) 


which is orders of magnitude more stringent for and Ut than laboratory 
limits. The same condition also allows only certain ranges of lifetimes for 
unstable neutrinos [12]. Radiative decay channels are additionally restricted 
because the created photons would otherwise influence the thermodynamic 
evolution too strongly. Other decay channels might be possible. The most 
common ones discussed for heavy neutrino decays are 




- 


vh - 

l^L+1 

(128) 

/+/- 

1 = e, p 

(129) 

^ I^L 4 

-vl+ 1^1 

(130) 

VH - 

^ l^L + X 

(131) 


where y corresponds to a light scalar like the majoron. Bounds on the radiative 
decay mode exist from reactor experiments [291], LAMPF [292] and from SN 
1987A [293]. The decay mode involving e’''e“-pairs is restricted in the region 
1-8 MeV by reactor data [294] and at higher energies by accelerator searches 
(see [295]). 

Experimental detection of the neutrino background will be extremly difficult. 
One of the suggestions is to take advantage of coherence in z/A^-scattering 
and by using cryogenic detectors, but also this proposal is far from practical 
realisation. A chance might be that eq.l25 is modihed during decoupling in the 
early universe by incomplete annihilation and hnite temperature QED-effects, 
which might increase by about 1 % [296]. This effect might be detectable 
in the future satellite missions MAP and PLANCK by its influence on the 
anisotropies in the cosmic microwave background [297]. 


8 Conclusions and outlook 


The question whether neutrinos have a non-vanishing rest mass influences 
research areas from particle physics up to cosmology, but is still an open 
question. At present all hints for neutrino masses are connected with neutrino 
oscillation effects namely the solar neutrino dehcit, the atmospheric neutrino 
anomaly and the evidence from LSND . The involved Am^ scales are in the 
MSW-(vacuum)solution, and eV-region respectively. 

Several theoretical models have been developed to describe these evidences in a 
consistent way [298,299]. With the appearance of the new CHOOZ and Super- 
Kamiokande results most of these models seem to be ruled out. Moreover 
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to explain all data in addition to the three standard neutrinos a new sterile 
neutrino seems necessary [300-303]. The scheme of Barger et ah [302] proposes 
a nearly degenerate and Vr in the eV-range and much lighter and z/5. 
The splitting between the eV-states and the light ones is determined by LSND 
data. The atmospheric anomaly can be solved by - Vr oscillations and both 
can act as hot dark matter. The solar neutrino problem is solved by Ve - 
oscillations which require the vg to be slightly heavier than z/g . 

The right answer will hopefully be given in the future, because of an increasing 
number of new experiments. 
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